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ABSTRACT 

Objects with radii of 10 pc to 100 pc and masses in the range from 10 6 M to 10 8 M Q 

i 1 have been discovered during the past decade. These so-called ultra compact dwarf 

^ ■ galaxies (UCDs) constitute a transition between classical star clusters and elliptical 

Q-C galaxies in terms of radii, relaxation times and U-band mass-to-light ratios. Using 

q , new data, the increase of typical radii with mass for compact objects more massive 

!_i ■ than 10 6 M Q can be confirmed. There is a continuous transition to the typical, mass- 

' independent radii of globular clusters (GCs). It can be concluded from the different 

relations between mass and radius of GCs and UCDs that at least their evolution must 
have proceeded differently, while the continuous transition could indicate a common 
formation scenario. The strong increase of the characteristic radii also implies a strong 
^ , increase of the median two-body relaxation time, t re \, which becomes longer than a 

£f-^ • Hubble time, th, in the mass interval between 10 6 Mq and 10 7 M Q . This is also the 

mass interval where the highest stellar densities are reached. The mass-to-light ratios 
J*"-. ■ of UCDs are clearly higher than the ones of GCs, and the departure from mass-to- 

<^> ' light ratios typical for GCs happens again at a mass of w 10 6 M Q . Dwarf spheroidal 

. galaxies turn out to be total outliers compared to all other dynamically hot stellar 

■ systems regarding their dynamical mass-to-light ratios. Stellar population models were 

' consulted in order to compare the mass-to-light ratios of the UCDs with theoretical 

. predictions for dynamically unevolved simple stellar populations (SSPs), which are 

!LJ ' probably a good approximation to the actual stellar populations in the UCDs. The SSP 

. £h \ models also allow to account for the effects of metallicity on the mass-to-light ratio. It 

. is found that the UCDs, if taken as a sample, have a tendency to higher mass-to-light 

$_i ' ratios than it would be expected from the SSP models assuming that the initial stellar 

mass function in the UCDs is the same as in resolved stellar populations. This can be 
interpreted in several ways: As a failure of state-of-the-art stellar evolution and stellar 
population modelling, as a presence of dark matter in UCDs or as stellar populations 
which formed with initial stellar mass functions different to the canonical one for 
resolved populations. But it is noteworthy that evidence for dark matter emerges only 
in systems with t re ] > th- 
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1 INTRODUCTION 

Star clusters can be defined as stellar population with a me- 
dian two-body relaxation time, t re i, shorter than a Hubbl e 
time, th, while galaxies would have t Te l > th l|Kroupall 19981 ). 
The dynamical evolution of the former is well described by 
pure Newtonian dynamics, while for the successful repre- 
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sentation of the latter either a significant amount of dark 
matter (DM) is required for Newtonian gravity to remain 
valid, or modified gravity needs to be invoked. By moving 
from two-body relaxation dominated systems to such where 
two-body relaxation plays no role, we thus observe the ap- 
pearance of fundamentally new physics. A transition class 
of objects between classical star clusters and galaxies may 
shed insights to the possible nature of the deviant dynamics 
apparent on galaxy scales. 

It has been almost 10 years since iHilker et all (|l999h 
and lDrinkwater et al.l (|2000l ) discovered these transition ob- 
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jects in the Fornax galaxy cluster. With apparent U-band 
magnitudes of < 19.5 mag at that distance, they can in 
principle be detected without difficulty. However, they can- 
not be discriminated from point sources with ground-based 
telescopes, except with the ones with the highest currently 
available resolutions. Because of this combination of small 
extension and high brightness they were usually thought to 
be foreground stars. Only a radial velocity survey of all ob- 
jects with a certain brightness in an area around the cen- 
tral galaxy of the Fornax cluster w as able to revea l their 
membership to that galaxy cluster. IPhillipps^ et all (120011 ) 
were the first ones to call them ultra compact dwarf galaxies 
(UCDs) , a ter m which is widely in use fo r this type of objects 
at the present. iDrinkwater et all (|2003t ) reported that these 
objects are not only distinct from the globular clusters in the 
Milky Way (MWGCs) by their higher U-band (L v ) lumi- 
nosity, but also by their larger radii and higher dynamical V- 
band mass-to-light (M/Lv) ratios. At the same time, there 
is no g ap in luminosity between g lobular clusters (GCs) and 
UCDs l|Mieske et al.ll2002l.l2004l ). lHasegan et al. I (|2005h dis- 
covered in the Virgo cluster massive compact star clusters 
with similar properties like the ones in the Fornax cluster, 
but called them dwarf-globul ar transition objects (DGTOs). 
Like lDrinkwater et al.1 (|2003l ). they state that the dynamical 
M/Lv ratios of some of the objects they disco vered are sig- 
nifican tly higher than the ones of the MWGCs. iMieske et all 
( 2006) concluded from the H/3 indices of UCDs in the Fornax 
cluster that they a re mo st likely of intermediate age, while 
lEvstigneeva et all ([2007) found the H/3 indices of UCDs in 
the Virgo cluster most consistent with old ages. Their stel- 
lar population has evolved passively for a long time in any 
case, which makes UCDs similar to most GCs and elliptical 
galaxies in this respect. 

Several formation scenarios that account for the physi- 
cal properties of the UCDs have been proposed: 

(i) UCDs are the mergers of many massive young clusters 
that formed in a star burst triggered by a galaxy-galaxy en- 
counter (e.g. like in the Antennae). After ~ lOGyr of dy- 
namical (an d stellar) evolution, such an object would resem- 
ble a UCD (|Kroupa(ll99g| ; iFellhauer fc Kroupall2002i). 



ii) UCDs are the most luminous GCs ( Mieske et al.l 



120021 ) 

(iii) UCDs are the central parts of nucleated galaxies that 
were disrupted by tidal forces as they moved in the grav- 
itational field of a larger galaxy. O nly the tightly bound 
cores survived until the present times (IZinnecker et al.ll 19881 ; 
iBassino et allll994 iBekki et al.ll2003l ; iGoerdt et al-lbOOsT T 

(iv) UCDs are the remna nts of the fundamental building 
blocks in galaxy formation l|Drinkwater et al"1l2004l ). 

Some bright UCDs in the F ornax cluste r and the 
Virgo cluster have been a nalysed bv lHilker et al. I (120071 ) and 
lEvstigneeva et ah! (|2007l ) very recently. They provide de- 
tailed high-quality data for 11 UCDs with dynamical masses 
between 10 7 Ma and 10 s Mq. Similar data have been ob- 
tained bv lReikuba et all (|2007l ) for compact objects in Cen- 
taurus A, but mostly with mass es between 10 6 Mq a nd 
10 7 Mq . They enlarge a sample by lHasegan et al.1 (|2005l ) in 
the Virgo cluster by 20 objects in the same mass range. 
Taken together, these data allow us to analyse the change 
of the internal parameters of massive comp act objects with 
mass or luminosity in more detail than IDrinkwater et all 



(2003) or lHasegan et all |2005l ). Furthermore, a comparison 
to other dynamically hot stellar systems (i.e. stellar systems 
whose stars are on randomised orbits) becomes possible, 
since samples with similar measured quantities are avail- 
able as well. The quantities that are considered here include 
their M/Lv ratio, Ty, and their projected (effective) half- 
light radius, r e , in dependency of their dynamical mass. 

Especially the dynamical M/Lv ratios of the 
UCDs has caugh t the attention of astronomers lately. 
lEvstigneeva et all (|2007l ) find the UCDs in their sample to 
be consistent with predictions fro m simple stel l ar pop ulation 
(SSP) models within the errors. iHilker et ail (|2007l ) note a 
tendency of the SSP models to under-predict the M/Lv ra- 
tios if a stellar population consist ent with observations i n the 
solar neighbourhood is assumed. lHasegan et al.1 j|2005l) find 
that some of the stellar systems they discuss have M/Lv ra- 
tios that imply extreme stellar populations in these objects. 
They suggest a presence of DM in these objects, provided 
that they are in dynamical equilibrium. This contradicts sce- 
nario (i), in which UCDs form DM free. Also if UCDs are 
nothing but very luminous GCs (scenario ii), they would be 
expected to be DM free, since G Cs of usual size are. The 
simulations by IBekki et al.1 (§003) on scenario (iii) predict 
DM free UCDs, since the DM halo of the progenitor galaxy 
of the UCD is found to be disrupted by the tidal interactions 
with the host galaxy of the UCD. This stands in con t rast t o 
the results from similar simulations bv lGoerdt et al.1 (2008), 
who found that a UCD can still be DM dominated if it is 
the stripped nucleus of a nucleated galaxy. Scenario (iv) also 
suggests dark matter in UCDs. A detailed analysis of the 
M/Lv ratios of the UCDs and their comparison to different 
SSP models may therefore give insights on their origin. 

The stellar population of the UCDs obviously plays a 
decisive role for the M/Lv ratio that has to be expected. 
The stellar population of each stellar system is determined, 
aside from an influence by stellar and dynamical evolution, 
by the stellar initial mass function (IMF), £(m), 



dN oc £(m) dm, 



(1) 



where m is the stellar initial mass and dN the number of 
stars in the mass interval [m,m + dm]. The IMF has to 
be distinguished from the present day stellar mass function 
(PDMF) which gives the number density of stars in depen- 
dency of stellar masses today. The IMF is a very useful con- 
cept, especially for a dynamically unevolved stellar system, 
because the number of stars that formed in the mass interval 
[m, m + dm] is conserved with time on the whole domain of 
the IMF. As a consequence, the PDMF and IMF are very 
similar for stars still on the main sequence at the present 
time. It turns out in Section 13.21 that UCDs can indeed be 
considered as dynamically unevolved stellar systems due to 
their mass and extension and therefore long relaxation time. 

In the past, there have been numerous efforts to in- 
fer the shape of the IMF from the PDMF as observed 
in resolvable stellar populations. There is common agree- 
ment that these observa tions are compa tible with the IMF 
originally proposed by ISalpeterl (Il955l ) for field stars in 
the solar neighbourhood: £(m) tx m~ a with a = 2.35 for 
0.4 Mq < m < 10 Mq. Later observations indicated that a 
is constant up to the highest obse rved stellar masses (which 
are between 120 Mq and 200 Mq , IWeidner fc Kroupal 12061 : 



lOev fc Clarke] 12003 : iFigerl 120051 ). but gets smaller below 
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0.5 Mq (|Kroupall2o"oil and references therein). The IMFs we 
consider for the stellar populations of the UCDs are guided 
by these results. 

With masses between 10 7 M and 10 s M and half-light 
radii mostly below 50 pc, UCDs may have formed contain- 
ing, within no more than some ten pc, between 10 s and 
10 6 O-stars or an order of magnitude more if the IMF was 
top-heavy. This is a scale of star formation beyond current 
theoretical reach, and it is therefore interesting to study the 
stellar content of these objects to probe the very extreme 
physics of their formation. 

Let us stress the importance of dynamical mass esti- 
mates for a meaningful discussion of the M/Ly ratios. This 
puts a hard constraint on the UCDs that can be included in 
this discussion since it requires high-resolution spectroscopy 
of faint objects. However, a dynamical mass estimate is in- 
dependent from the total luminosity of the stellar system. 
Instead, the mass estimate is based on the surface brightness 
profile and the wid t h of t he spectral lines as described in de- 
tail in Iffilker et al l (|2007l ). Dynamical mass estimates clearly 
rely on a number of assumptions that cannot be verified 
easily, but mass estimates for unresolved stellar populations 
based on stellar population models do so as well. The true 
advantage of the dynamical mass estimates for this work is 
that they allow an independent estimate for the M/Lv ratio 
that can be compared to theoretical predictions from stellar 
population models. 

This paper is organised as follows. In Section[2]a sample 
of different dynamically hot stellar systems, including UCDs, 
is introduced. Section [3] is dedicated to the dependencies of 
internal parameters of dynamically hot stellar systems on 
their mass. The M/Lv ratio of UCDs, GCs and elliptical 
galaxies is compared to the predictions from simple stellar 
population models in Section [4] While doing this, we take 
the influence of their metallicity on their luminosity into ac- 
count. Section [5] contains a discussion of the transition from 
GCs to UCDs. Furthermore the reliability of our results con- 
cerning the M/Ly ratio of UCDs is addressed. We conclude 
with Section [6] 



2 THE DATA 

One of the tasks performed in this paper is to compare UCDs 
to other dynamically hot stellar systems. This requires a 
set of data which spans over many orders of magnitude in 
dynamical mass. A homogeneous data sample is unfortu- 
nately not available due to the diversity of the objects. We 
therefore collect data from different sources in the litera- 
ture, where comparable parameters have been measured or 
where at least a correlation between the measured data to 
the ones that are to be compared is known. In the following, 
we specify the sources for our data and how we derived the 
quantities we use in this paper from them, if necessary. 

2.1 Massive Compact Objects 

It is convenient in this paper to introduce massive compact 
objects (MCOs) as a collective term for all stellar systems 
in the sample discussed here that should neither be denom- 
inated as MWGCs nor as elliptical galaxies. This definition 
of MCOs includes a number of objects that are considered 



as UCDs in other works. The motivation for the introduc- 
tion of this term lies in the fact that the sample of objects 
discussed here also includes a number of objects which in 
their entirety seem to mark a transition between GCs and 
UCDs. This will become apparent below. A clear distinction 
between GCs and UCDs is thereby problematic here. 

We differentiate the MCOs by the way their dynamical 
masses were estimated: 

For the 19 MCOs listed in Tab. [T] the mass estimate 
included the fitting of a density profile to each one of them 
individually. These 19 objects are 12 MCOs from the Virgo 
cluster, five UCDs from the Fornax cluster as well as two 
objects from the Local Group: u> Cen in the Milky Way and 
Gl in Andromeda. We consider uj Cen as an MCO instead 
of an MWGC because of its spread in [Fe/H], which sets it 
apar t from every other s t ar cluster in the halo of the Galaxy 
(e.g. iKavser et aD l200d : IVillanova et al.1 l2007t ) We refer to 
them as "MCOs with mass distribution modelling" . 

We also i nclud e 20 objects in Centaurus A from 
iReikuba et ail (|2007l ) for which measurements of the veloc- 
ity dispersion and at least one colour index are available. 
Tab.[2]lists their properties. Their mass i n Ma is calcul ated 
by using a virial mass estimator given in ISpitzerl fll987h : 

M a ~ 10G~V c cr 2 , (2) 

where r e is the projected half-light radiuij]] in pc and a is the 
global velocity dispersion in pc Myr" 1 . G is the gravitational 
constant, which is 0.0045 pc 3 Mg 1 Myr -2 We refer to them 
as "MCOs with global mass estimate". 



2.2 Globular clusters 

We compare the MCOs to the MWGCs for which 
McLaughlin & van der Marell (|2005h calculated dynamical 
M/Lv ratios (listed in their table 13). Their value for the 
effective half-mass radius and their estimate of the dynami- 
cal M/Lv ratio in the U-band for the King model is used in 
this wor k. By u s ing th e absolute magnitude in the U-band 
given in Irlarria (| 1996h . the cluster mass can be calculated 
from its M/Lv ratio. 

It can hardly be expected that such a limited sample is 
representative for GCs in general. Nevertheless, this seems 
to be the case to some extent, as surveys of extragalactic 
GC s yst ems show (e.g. lLarsen et al.l 1200 ll : IChandar et alJ 
2004 and Jordan et al. 2005 concerning the radii of GCs, and 
iRicht lcr 2003 and Jordan et al.l 120071 concerning the absolute 
magnitudes of GCs, which indicate their masses if a constant 
M/L ratio for them is assumed) . It therefore seems possible 
to take the distribution of the radii and the masses of the 
MWGS as a rough representation of GCs in general. The 
advantage of the chosen sample is that, as for the MCOs, 
mass estimates from velocity dispersions are available for 
them. 



1 Actually, it is the half-mass radius that enters into eq. {2J , but 
we assume that the mass density follows the luminosity density 
whenever necessary. This allows us to identify the half-mass ra- 
dius with the half-light radius. 
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Table 1. Properties of MCOs with masses from mass distribution modelling. The contents of the columns are the following. Column 1: 
The name given to the MCO (the same as in the source papers), Column 2: The projected half-light radius of the MCO, Column 3: The 
global velocity dispersion of the MCO, Column 4: The central velocity dispersion of the MCO, Column 5: The absolute magnitude of the 
MCO in the V-band, Column 6: The dyn amical mass of the MCO, Col umn 7: The M/Ly ra tio of the MCO, Colum n 8: Refer ences to the 
paper s th at are the basis fo r our d ata: 1: lEvstignccv a et al. (2007), 2 : lHasegan et al. 3: Hilk er et al. I ll2007h , 4: iBaumgardt et al,l 

112003! ), 5: Ivan de Ven et ahl 11200611 , 6: lHarrisI dl996Tl . Some errors are marked with an asterisk; they have not been published so far. 



Name 


r e 




a 




My 


M 




M/L v 


Rcf 




[pc] 




[kixi S J 


[Km S J 


[mag] 


[1O 6 M ] 


[IVlQ 1j q J 




VUCD1 


11.3 ± 


0.7* 


32.2 ± 2.4 


39.3 ± 2.0 


-12.26 


28.0 ± 


5.0 


4.0 ±0.7 


1 


VUCD3 


18.7 ± 


1.8* 


35.8 ±1.5 


52.2 ± 2.5 


-12.58 


50.0 ± 


7.0 


5.4 ±0.9 


1 


VUCD4 


22.0 ± 


2.7* 


21.3 ±2.0 


26.9 ± 2.3 


-12.30 


24.0 ± 


6.0 


3.4 ±0.9 


1 


VUCD5 


17.9 ± 


0.8* 


26.4 ±1.6 


32.5 ± 2.3 


-12.32 


29.0 ± 


1.0 


3.9 ±0.6 


1 


VUCD6 


14.8 ± 


3.1* 


22.3 ±1.8 


29.6 ± 2.2 


-12.10 


18. 0± 


5.0 


2.9 ±0.9 


1 


VUCD7 


96.8 ±20* 


27.2 ±4.6 


45.1 ± 1.5 


-13.44 


88.0 ± 


21.0 


4.3 ±1.1 


1 


S417 


14.36 ± 


0.36 


26.4 ±2.7 


31.7 ± 1.4 


-11.78 ±0.16 


27.0 ± 


5.0 


6.6 ±1.5 


1.2 


UCD1 


22.4 ± 


1.0 


27.1 ± 1.8 


41.3 ± 1.0 


-12.19 


32.1 ± 


3.6 


4.99 ±0.60 


3.1 


UCD2 


23.2 ± 


1.0 


21.6 ±1.8 


31.3 ± 0.6 


-12.27 


21.8 ± 


3.1 


3.15 ±0.49 


3.1 


UCD3 


89.9 ± 


6.0 


25.0 ±3.4 


29.3 ± 1.2 


-13.57 


94.5 ± 


22.0 


4.13 ±0.98 


3.1 


UCD4 


29.6 ± 


2.0 


22.8 ±3.1 


37.3 ± 0.6 


-12.45 


37.3 ± 


8.6 


4.57±1.11 


3.1 


UCD5 


30.0 ± 


2.5 


18.7±3.2 


28.7 ±0.8 


-11.99 


18. 0± 


5.0 


3.37±0.85 


3.1 


S314 


3.23 ± 


0.19 




35.3 ± 1.4 


-10.91 ±0.16 


5.8± 


1.0 


2.94 ±0.68 


2 


S490 


3.64 ± 


0.36 




42.5 ± 2.7 


-11.00 ±0.16 


8.7± 


2.1 


4.06 ±1.15 


2 


S928 


23.16 ± 


1.37 




22.4 ± 1.0 


-11.58 ±0.16 


21.3 ± 


2.9 


6.06 ±1.23 


2 


S999 


20.13 ± 


0.98 




25.6 ± 1.4 


-11.08 ±0.16 


21.6 ± 


2.9 


9.36 ±1.87 


2 


H8005 


28.69 ± 


0.55 




10.8 ± 2.3 


-10.83 ±0.16 


5.5± 


2.3 


2.98 ±1.35 


2 


Gl 


8.21 






25.1 ± 1.7 


-10.94 


8.2 ± 


0.85 


4.10 ±0.42 


4 


ui Cen 


6.70 ± 


0.28 


16.0 


19.0 ± 1.5 


-10.29 


2.5 ± 


0.1 


2.5 ±0.3 


5.6 



Table 2. Properties of the compact object s in Cen t aurus A. Here the mass was calculated by using the same mass estimator for all objects, 
namely eq. J2j. All d ata are from iReikuba et al.l J2007tl . The contents of the columns are the following. Column 1: The identification 
of the object (like in lReikuba et al .1 l|2007l l'). Column 2: The effective (projected half light) radius of the MCO, Column 3: The global 
velocity dispersion, Column 4: The estimated (dynamical) mass, Column 5: The M/Ly ratio. 



Name r c a M„ M/Ly 

[pc] [krns" 1 ] [1O 6 M ] [MqL" 1 ] 



HGHH92-C7 


7 


5 ± 


1 


21 


fi +1.0 

°-2.6 


7 


O+0.7 
°-1.9 


3 


3+? 


HGHH92-C11 


7 


8±0 


1 


19 


g+0.9 
°-2.3 


6 


7 +0.6 
'-1.6 


5 


7±i 


HHH86-C15 


5 


3±0 


7 


11 


1+0.7 

x -0.7 


1 


c+0.2 
°-0.5 


2 


q + 


HGHH92-C17 


5 


7± 


1 


20 


Q+1.6 

a -1.6 


5 


c+0.5 
8 -1.4 


3 


8+; 


HGHH92-C21 


7 


0±0 


1 


19 


q+0. 8 
J -2.3 


5 


c+0.5 
8 -1.4 


1 




HGHH92-C22 


3 


8±0 


1 


17 


Q+0.1 

a -0.1 


2 


O+0.3 
8 -0.7 


3 


o+; 


HGHH92-C23 


3 


3±0 


1 


31 


q+1.4 
^-3.9 


7 


9+0.7 
-1.8 


1 




HGHH92-C29 


6 


9±0 


1 


16 


1+0.8 
1 -0.8 


1 


1+0.4 
^-1.0 


4 


4±i 


HGHH92-C36 


3 


6±0 


3 


15 


7+1-9 
' -1.9 


2 


+0.3 
U -0.6 


2 


g+0 
D -0 


HGHH92-C37 


2 


9±0 


3 


12 


g+0.8 
°-0.8 


1 


i+O.l 
-0.3 


1 


5 +0 
°-0 


HHH86-C38 


2 


8±0 


2 


14 


n+l.l 
z -l.l 


1 


q+0.2 
^-O^ 


1 


O + 

8 -0 


HGHH92-C41 


4 


5 ± 


1 


11 


c+1.3 
°-1.3 


1 


4+0-1 
^-0.3 


2 


9+0 


HGHH92-C44 


5 


7±0 


1 


13 


i+l.O 
-1.0 


2 


q+0.2 


3 


9+° 


HCH99-2 


11 


4± 1 


1 


14 


1+0.5 
L -0.S 


5 


q+0.7 
°-1.5 


4 




HCH99-15 


5 


9±0 


2 


21 


o+1.7 
°— 1.7 


6 


9+0.6 
z -1.5 


3 


4 +0 


HCH99-16 


12 


1±0 


6 


9 


c + 1.4 
°-1.4 


2 


c+0.3 
°-0.6 


2 




HCH99-18 


13 


7±0 


3 


21 


9+11 
Z -l.l 


14 


q+1.3 
°-3.5 


1 


rtl 


HCH99-21 


7 


1±2 


7 


10 


g+2.3 
°-2.3 


1 


q+0.7 
M -1.0 


1 


7+0 
'-1 


R223 


2 


6±0 


3 


14 


4 +l.5 


1 


q+0.2 

°-0.4 


2 


q+0 


R261 


1 


9±0 


4 


14 


g+0.7 
"-0.7 


1 


+0.2 
U -0.3 


1 


1+° 

i -0 
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2.3 Early- type galaxies 

We also compare the MCOs to more massive dynamically 
hot stella r systems by making use of some of the data pub- 
lished bv lBender et al.l (|l992l ). i.e. their values for the cen- 
tral velocity dispersion, oo, the projected half-light radius, 
r e and the absolute magnitude in the B-band of elliptical 
gal axies and bulges of ea rly-type spiral galaxies in their sam- 
ple . IBender et all (|l992h give a simple formula for estimating 
the King mass from r c and an, which we use as well for the 
objects from their paper: 



5G 1 r e ao, 



(3) 

in pcMyr -1 and G = 



with r c in pc, <tq 
0.0045 p^M" 1 Myr" 2 . 

If these objects are to be compared to the MCOs, their 
V^-band luminosities have to be estimated from their S-band 
luminosities, since for the MCOs luminosities in the V-band 
are measured. It is known that there is a correlation be- 
tween the luminosity and the colour of elliptical galaxies. 
However, given the weakness of this dependency, we think 
that accounting fo r it (e.g. with the d ata on colour of the 
same galaxies from IBender et al"1 [T993) would probably not 
pay the effort. This becomes evident, if the uncertainties 
connected to the mass estimates from eq. © especially are 
considered (see Section 12.411 . Therefore, adopting a uniform 
B — V colour index of 0.9 seems a reasonable approximation 
for the purpose of this paper. 

To enhance t he sample, data on nucleated dwarf ellip- 
tical galaxies from lGeha et al.l l|2003h are included. 

Data on dwarf spheroida l galaxies (dSphs ) are a lso in- 
cluded. They are taken from iMetz fc Kroupal (|2007l ). their 
table 2, because their data on dSph s are more up to date 
than the ones in IBender et all <|1992T >. The half-light radii of 
the dSphs are not listed in that table, but are usually found 
in the references given there (with the exception of And II, 
for which the half-light rad ius is taken from the paper by 
iMcConnachie fc Irwirj|2006h . 



2.4 A note on different dynamical mass estimators 

The dynamical mass of each of the objects introduced above 
was estimated in one of three different ways. While for some 
objects the mass estimate included the fitting of an individ- 
ual density profile to them, for other objects the mass was 
calculated by using one of two global mass estimators. The 
choice of the mass estimator depended on whether a or oq 
of the a stellar system was measured. This raises the ques- 
tion whether the mass estimates obtained in these different 
ways are indeed comparable. If they are comparable, two 
requirements should be fulfilled: 

(i) There should not be a tendency for one method to 
over- or underestimate the mass. 

(ii) Applying different mass estimators on the same ob- 
ject should give similar results. 

This can be tested on the 19 MCOs in Tab. [[] where 
a and on or ao only is available beside the mass estimate 
using an individual density profile, M, which is probably the 
most reliable one and therefore is considered as a standard 
here. Fig. [1] shows the masses as determined by using the 
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Figure 1. Plot of the ratios between global estimates and mass 
estimates including mass distribution modelling for the 19 MCOs 
in Tab. [T] against the estimate for their mass from an individual 
fit. Open squares show M&/M and circles show M^q/M. 



global mass estimators in comparison to the mass from an 
individual density profile fit. 

As a measure for the mean deviation of the mass es- 
timated using eq. ((2j), M CT , and the mass estimated using 
eq. Q, Mao, from M, we calculate AM, = y"]^ 1 



I Mi 



M ai \ and AM CT 



where Ni and 



N2 denote the number of objects that are included for that 
summation. This results in AM„ = 8.5 x 10 6 Mq for the av- 
erage deviation of M a from M. This value can be compared 
to the mean value for the mass M of the same MCOs, with 
the masses as they are estimated using individual models 
for the density profile, which is M = 36.2 x 1O 6 M . This 
means that the average deviation of M a from M is about 
23%. 

Similarly, the average deviation of M CT o from M can be 
calculated: AM a o = 12.5 x 1O 6 M0. If this is again compared 
to M of the according MCOs, it turns out that the average 
deviation of M a o from M is about 44%. The larger discrep- 
ancies between M and M a o than between M and M a is at 
least partially due to the uncertainties to the inner density 
profiles of the MCOs, because the central structure of an 
MCO strongly influences the value that is determined for 
its on. 

The (relative and absolute) discrepancy between M and 
M a or M ct0 is the largest for VUCD7. However, VUCD7 is 
one of those MCOs that are best fit by a two-component 
(King+Sersic) density profile, in contrast to most of the 
other MCOs. It is therefore not surprising that the mass 
estimators eq. ([2]) and eq. Q fail here, since they assume a 
King profile. This illustrates the risk connected to assuming 
a single typical profile for a number of objects. Excluding 
VUCD7, the average deviation of M a from M can be low- 
ered to about 10%, and the average deviation of M^o from 
M can be lowered to about 24%. 

In summary, the three ways to estimate the dynam- 
ical mass seem to produce comparable results. Note that 



6 Dabringhausen et al. 



also iHilker et all (|2007l ) and lEvstigneeva et all <|2007f > usu- 
ally find that the internal parameters derived from global 
King estimators (a = 2) are almost identical to the pa- 
rameters derived using mass distribution modelling. We will 
therefore not discriminate between M a , M a o and M any 
further, but denote all dynamical masses as M. 



3 DEPENDENCIES ON DYNAMICAL MASS 

In this section, the effective radii, median relaxation times, 
central densities and M/Ly ratios are compared to each 
other. 



3.1 Dependency of the effective radius on mass 

In Fig. [2] the mass dependency of r c of the MCOs and other 
dynamically hot stellar systems is plotted. Some well estab- 
lished observations can be identified easily in this plot: The 
strong correlation b etween M and r e for elliptical galaxies 
l|Bender et al.|[l992i ) in the high mass ran ge and the absence 
of a dependency o f r c on M for GCs (IMcLaughlinl |2000| ; 
IJordan et ail koOS) at the lowest masses. Remarkable is the 
large spread of radii at intermediate masses which becomes 
largest in the mass interval of 10 7 M Q < M < 10 8 M Q , the 
mass interval where the rather compact UCDs as well as 
the (typically about an order of magnitude) more extended 
dSphs lie. The underlying assumption for this statement is 
that dSphs are objects in (or close to) vir ial equilib rium. 
This has been argue d to be the case by e.g. IWul (|2007i ) and 
iGilmore et al.l (|2007t ) for at least those dSphs that are most 
distant to the Galactic centre, although thi s would imply 
extrem ely high M/Lv ratios in some cases. IGilmore et al.l 
state that there is a bimodality of the characteristic 
radii of objects in the mass range 10 7 Mq < M < 10 8 Mq, 
i.e. an almost complete absence of objects with r c ~ 100 pc. 
In Fig. 21 they are indeed only represented by VUCD7 and 
UCD3j (and M32 at a higher mass). One way to interpret 
this is to consider UCDs and the dSphs as two kinds of 
stellar system s that formed under different conditions, as 
IGilmore et al. I (120071) propose- 
However. iMetz fc Kro upa (2007]) argue that the forma- 
tion of dSphs may have been triggered by the tidal forces 
in an encounter between two galaxies, i.e. they propose in 
principle the same sc e nario for the formation of dSphs which 
iFellhauer fc Kroupal (|2002l ) suggested for the formation of 
UCDs. The morphological differences can be understood 
in terms of the influence of the surroundings on the star- 
forming regions: the dSphs can form from star cluster com- 
plexes in a weak tidal field (e.g. the tidal arm of the Tadpole 
galaxy), while the UCDs form in a strong tidal field (e.g. 
the Antennae galaxy). This scenario is supported with the 
observation that the orbital angular momenta of the satel- 
lite g alaxies of the Milky Way are correlated (jMetz et alj 
120081 ). It can also offer an explanation for the seemingly 



high M/Lv ratios of some of the dSphs, if they are largely 



unbound phase-space structures and therefore cann ot be de- 
scrib ed by simple application of Jeans' equations (jKroupal 
ll997l V 

It is surprising that the MCOs lie on the same rela- 
tion between mass and radius as massive elliptical galax- 
ies with masses > 10 Mq, while elliptical galaxies with 
lower masses (i.e. objects in the intermediate mass range) 
mostly lie on a different relation, which points towards 
the parameter space of dSphs. This could be evidence for 
the low-mass ellipti cal galaxies being mo s tly of tidal ori- 
gin, as p roposed by Okazaki fc TaniguchU (|2000l ) (also see 
fig. 7 in iMonreal-Ibero et al.l 120071 ). and as discussed by 
IMetz fc Kroupal 120071 for dSphs. The few compact low- 
mass elliptical galaxies can then be interpreted as low-mass 
counterparts of the elliptical galaxies more massive than 
> 10 11 Mq . 

Following the above interpretation, some objects are 
thus excluded for quantifying the relation between mass and 
radius that MCOs share with massive elliptical galaxies in 
a least squares fit. These objects are, b esides the MWGCs 
and the dSphs, the d warf ellipticals fromlGeha et al.l (120031 ) 
and the galaxies that lBender et al.l (| 1992T) define as "bright 
dwarf ellipticals'ff]. The exclusion of the latter two groups 
may seem somewhat arbitrary, but it turns out that they de- 
fine the apparent turn-off from the relation for the remaining 
objects (i.e. bright elliptical galaxies, galaxy bulges, compact 
ellipticals and MCOs) at > 10 11 Mq quite well. Assuming a 
function of the form 



M 



(4) 



2 Note that IGilmore et al. I 1120071) c onside r a half-light radius of 
only 22 pc (from IpTmkwater et alj |2003)) for UCD3 and omit 
VUCD7 from their discussion. 



PC a \ 10 e Mq 

for the relation between M and r e , which corresponds to a 
straight line in Fig. [2] leads to 

a = 2.95^Q22> 
b = 0.596 ±0.007, 

for the best-fitting parameters. If the MCOs are not used 
for the fit, 

O r ,,+0.91 

b = 0.608 ± 0.025, 

is obtained, i.e. within the errors the same relation as with 
the MCOs. The small impact that excluding the MCOs has 
on the fit is demonstrated in Fig. [5] by plotting eq. @ with 
both sets of values for a and b. This veryfies that the MCOs 
lie along the same relation between M and r c as massive 
elliptical galaxies. 

For comparison, an analogous fit to all elliptical galaxies 
as well as the dSphs (but without the MCOs) is performed. 
This corresponds to the hypothesis that these objects are 
drawn from a homogeneous population, which obeys a single 
relation between mass and radius. This leads to 

a = 34.8±S;5, 

b = 0.399 ±0.019, 

for the best-fitting parameters. However, the distribution of 
the massive elliptical galaxies is clearly asymmetric around 



3 i.e. those galaxies which have My > —18.5 and are no t classi- 
fied as "compact dwarf ellipticals" bv lBender et al.l (Il992h 
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Figure 2. Plot of the half-light radius, r c , against mass, M, for different types of dynamically hot stellar systems. The symbols that are 
used have the following meaning: Open circles for MWGCs, open diamonds for MCOs with global mass estimate (i.e. calculated from 
eq. f5J, fille d circles for MCOs with (the probably mor e reliable) mass estimate s from mass distribution modelling (i.e. mass estimates 
taken from lHasegan et~alll2005l. iHilker et al]|2007l and lEvstieneeva et al.ll2007l as well as the mass estimates for w Cen and Gl), open 
squares for dSphs, triangles for elliptical galaxies and filled diamonds for bulges of early- type spiral galaxies. Errors are comparable to 
the symbol sizes. The lines show fits to the data for a relation between mass and radius for bright ellipticals, compact ellipticals and 
bulges (dashed line), bright ellipticals, compact ellipticals, bulges and MCOs (solid line) and all elliptical galaxies, bulges and dSphs 
(dashed-dotted line). Most elliptical galaxies with low brightness have been excluded from the first two fits, see text for more details. 
They are marked with a cross. Note that the underlying assumption for the mass estim ates is that the stellar systems are essentially 
undisturbed by tidal fields, which may be wrong for the dSphs especially (Kroupa 1997). 



this relation, which suggests that the first two relations are 
a better fit to them. 

We note th at the larger sample of elliptical galaxies 
which is used bv lGraham et all l|200fj ) shows a very similar 
distribution of chara cteristic radii against mass, although 
Gra ham et al. I (|2006l ) estimated the masses of the galaxies 
differently to the approach chosen here, namely by assuming 
a stellar population for them and calculating their masses 
from their luminosities. 

The radii of MCOs are thus, unlike the ones of GCs, 
correlated to their masses. The comparison of the massive 
MCOs with the MWGCs shows that the characteristic radii 
of GCs are indeed typically about an order of magnitude 
smaller than the ones of the massive MCOs. However, Fig. [2] 
also seems to suggest a rather fluent transition between ob- 
jects that lie on the scaling relation for GCs and objects that 
lie on the scaling relation for elliptical galaxi es at a mass of 
about 10 6 Mq. This confirms the conclusions lHasegan et all 
(2005) have drawn based on fewer data. 

This change of typical radii cannot be due to an obser- 
vational bias against small radii for more massive objects, 



since MCOs are identified by their brightness, their mem- 
bership to a galaxy cluster and t heir compactness. The data 
on rather low-mass MCOs from lHasegan et~aH l|2005h and 
iReikuba et al.l (|2007T ) (both indicated as open diamonds in 
Fig- HI indeed include objects with radii on both scales. Con- 
sequently, this change of the typical r e must be connected to 
a difference in evolution or formation of objects less massive 
than « 10 6 Mq and more massive than ~ 10 7 M©. 



3.2 Dependency of the median two-body 
relaxation time on mass 

The median two-body relaxation time is closely connected 
with mass and characteristic radius of an o bject. It is given 
inMyr in a formula originally found by ISpitzer fc Hard 

CaziD, 



1 h 

GW 



(5) 



0.061 N 
tlcl ~ log(0.4JV) X 

where N is the number stars in the cluster, its half-mass 
radius in pc, M its mass in Mq and G the gravitational con- 
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Figure 3. The median relaxation time, t re i, plotted against dynamical mass, M. Contrary to Fig. [2] this figure shows MWGCs and 
MCOs only. The dashed line marks the current age of the universe. The symbols are as in Fig. [2] One MCO is plotted with typical errors. 



stant, which is 0.0045 pc 3 Mg 1 Myr -2 . t rcl can be considered 
as a measure for the relaxation time at the half mass-radius. 
However, eq. (f5]) is only an approximation: It is obtained un- 
der the assumption that the stars move in a smooth poten- 
tial and are only disturbed by two-body encounters (i.e. no 
binaries), beside the supposition that the cluster is in virial 
equilibrium. 

Eq. (0 includes parameters which are not known for 
most of the MCOs, but it can be transformed into one that 
only depends on M and the effective half-light radius, r c , as 
free parameters if some assumptions are made. It can then 
be applied to the data in this paper. This is done by assum- 
ing that the mass is dist ributed as th e luminosity and by 
substituting r c = 0.75 r h (|SpitzerJ | l9 87). We further assume 
a mean stellar mass of 0.4 M© in concordance with the mean 
stellar mass in a stellar population with the canonical IMF 
(see eq. 1 10p . This yields 



0.234 



log(M/M ) 



Mrj 
G 



(6) 



in the same units as eq. Jjj}. An inspection of eq. ((6)1 reveals 
that r c dominates the behaviour of t le \ due to its power. 
Therefore, a plot of t re i against M looks very similar to a 
plot of r c against M (Fig. [2}. 

In Fig. El trd is plotted against M of MWGCs and 
MCOs only. The stated similarity to Fig. [5] in the according 
mass range is apparent. The new and important piece of in- 
formation that can be read off Fig.[3]is how t re i of the objects 



compares to a Hubble time. It is clearly below a Hubble time 
for most MWGCs, while it is clearly above a Hubble time 
for all MCOs more massive than 10 7 Mq . This corresponds 
to the increase of the typical radii in the mass interval from 
1O 6 M to 10 7 M Q . As MWGCs and MCOs are considered 
to be old objects, this implies that MWGCs can have under- 
gone considerable dynamical evolution since their formation 
while massive MCOs have not. Consequently, massive MCOs 
are much less vulnerable to mass loss driven by two-body re- 
laxation. 



3.3 Dependency of the central density on mass 

It is worthwhile to consider the impact of the development 
of the typical radii with dynamical mass on the central den- 
sity of the MWGCs and MCOs. The central density is here 
defined as the mean density within the projected half-light 
(i.e. half-mass) radius. It is plotted in Fig. [4] against mass. 

The independence of the MWGC radii on their dynam- 
ical mass translates into an increase of the central density 
with dynamical mass. The increase of the typical radii above 
a dynamical mass of 10 6 Mq , as visible in Fig. [5] is strong 
enough for a slow decrease of t he central dens i ty to occur. 
It has already been noted by iBurstein et alj (|l997l ) that 
there is a maximum global luminosity density for early-type 
galaxies, which is proportional to M~ A ^ 3 . In this light, the 
decrease of the densities with mass for the MCOs is only a 
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Figure 4. The central density of MWGCs and MCOs plotted against dynamical mass. The symbols are as in Fig. \3\ The dashed lines 
indicate constant densities: assuming a mean stellar mass of 0.4 Mq, the lower dashed line indicates a density where the mean distance 
between stars is 6000 AU (about 100 times the diameter of the orbit of Neptune), and the upper dashed line indicates where the mean 
distance between stars is 3000 AU (about 50 times the diameter of the orbit of Neptune). 



consequence of the common relation between the MCOs and 
the massive elliptical galaxies that was found in Section |3. II 

3.4 Dependency of the M/Lv ratio on mass 

Fig. [5] shows the dynamical M/Lv ratios of the sample 
against the luminosity in the V-band (upper panel) and 
the dynamical mass (lower panel). It is visible from this 
figure that the dSphs with the lowest V-band luminosities 
al so have the h ighest M/Lv ratios, as was already noted 
in iMateol (| 19981 ). Other than that, the general distribution 
of the data in both panels is almost identical, except for a 
steeper rise of the M/Lv ratios from the MWGCs to the 
MCOs when they are plotted against Ly . 

It might be tempting to identify the gap in the lumi- 
nosity sequence at » 10 8 L© with the borderline between 
a star cluster-like population to the left and a galaxy-like 
population to the right. However, the homogeneous sample 
of faint early - type galaxies in the Fornax cluster observed by 
Hilkor et al.l (|200of ) does not show such a gap in luminosity 
down to the luminosities of dSphs. The gap visible in Fig. [5] 
is thus most likely an artefact caused by the inhomogeneity 
of our data sample. 

The spread of the M/Lv ratios of the dSphs is very 
striking in Fig. [5] With M / L v ratios of several 10 2 M /L , 
some of them are total outliers compared to all other dy- 
namically hot stellar systems. It is especially the spread of 



their M/Lv ratios that supports the notion that dSphs can- 
not be treated as objects in dynamical equilibrium. If they 
were in dynamical equilibrium, DM haloes with very differ- 
ent properties would have to be assumed for objects that are 
quite similar to each other as far as the properties of their 
baryonic matter are concerned. 

It can be seen for the remaining objects that almost 
every MCO above a mass of 10 6 Mq has a Mj Lv ratio which 
is manifestly higher than the mean value for MWGCs. As 
for the radii, the transition from the M/Lv ratios of GCs 
to the ones of MCOs seems fluent. The objects classified as 
some kind of elliptical galaxy (including bulges of early-type 
spiral galaxies) span the whole range of M/Lv ratios that is 
occupied by GCs and MCOs, with bulges and large elliptical 
galaxies having a larger spread to higher M/Lv ratios. 

It should be remembered in this context that the masses 
of the early-type galaxies that are used to determine their 
M/Ly-ratios have been calculated with eq. Q, i.e. the mass 
estimates are based on the distribution of the visible mat- 
ter. If these galaxies are embedded in DM haloes, the mass 
estimates are too low for the total masses of the galaxies, 
but are still good approximations for the mass of those parts 
of the galaxies that are dominated by baryonic matter. It is 
noteworthy that evidence for DM only emerges in objects 
with t IC \ > th (also see Fig. I10l and ll2[) . 

The physical reasons for the distribution of M/Lv ra- 
tios are a rather complicated issue. It mainly depends on 
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Figure 5. Dynamical M/Ly ratio plotted against luminosity in the V-band, Ly (upper panel), and mass, M (lower panel). The symbols 
are as in Fig. [2] The errors to the values of the MCOs are not much larger than the symbol size. 
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two things: A possible non-baryonic DM content in the ob- 
jects and the stellar populations of the objects. The M/Lv 
ratio of a stellar population is influenced by its star forma- 
tion history, its IMF, the metallicity of the stars and by how 
much the stellar population was altered by dynamical evolu- 
tion. Unfortunately, most of the objects in our sample can- 
not be resolved into stars so far, which makes it impossible 
to determine their stellar populations directly. Nevertheless, 
observations of these objects and theoretical considerations 
can give some clues on their stellar populations. Some of 
these findings are summarised below. 

• MWG Cs contain old ste l lar populations (older than 
» 10 Gyr, IVandenBerj |200d : ISalaris fc Weiss! l2002h . The 
MCOs in the Virgo c luster seem to have similar ages 
l|Evstigneeva et al.ll2007l ). bu t the MCOs in the F ornax clus- 
ter might be a bit younger (|Mieske et al.| [2006). The ages 
of elliptic al galaxies are found to range from a fe w Gyr to 
> 10 Gyr (|Trager et al.ll200d ; lAnnibali et al.ll2007l ). 

• MWGCs are known to have low metallicities. The 
metallicities of the MCOs are, if estimated, consistent with 
those of metal-rich MWGCs. Elliptical galaxies have about 
solar metallicities in their central parts I Trager et al.|[2000l ; 
lAnnibali et~al] 120071) and a decrease of their metallicities 
towa rds their outer regions (|Tantalo et al.l ll998; Bae s et ail 
12007ft . 

• Dynamical evolution can lower the M/Lv ratio of a 
stellar system noticeably, if the time scale for its dynamical 
evolution is shorter than the time scale for the evolution of 
its stars |Baumgardt fc MakindkoOSl ; iBorch et al.ll2007l ). 

With this information, the different M/Lv ratios of the 
objects plotted in Fig. [5] become understandable at least 
qualitatively. The rather low M/Lv ratio of MWGCs can be 
understood as an effect of their low metallicity and the con- 
siderable dynamical evolution that was sugg ested for them in 
sectio n [3721 Considering the lifetimes iBaumgardt fc Makinol 
(2003) expect for a sample of MWGCs (while accounting for 
the tidal field of the Galaxy) and their results for the devel- 
opment of the M/L ratio as a function of the star cluster 
lifetime, a decrease of the M/L ratio by about 0.1 Mq Lq 1 to 
0.3 Mq Lq 1 compared to the M/L ratio of a dynamically un- 
evolved stellar population would seem typical for MWGCs. 
The massive MCOs and the elliptical galaxies on the other 
hand are more metal-rich and due to their size and exten- 
sion dynamically almost unevolved. This might be able to 
explain higher M/Lv ratios compared to MWGCs even if 
they do not contain DM. Note however that a DM con- 
tent in elliptical galaxi es has been discussed: quite recently, 
ICappellari et al. 1 (120061 ) estimated a median DM content of 
~ 30% within the half-light radii of a sample of elliptical 
galaxies, if an IMF as in the Solar neighbourhood is as- 
sumecfl. The large spread of the M/Lv ratios of ellipticals 
is not surprising in the light of their large age spread. Also 



4 ICappellari et a l. (2006) do not discuss gas as a possible contrib- 
ut or to the non-luminou s matter. However, considering the results 
bv lCombes et aT I ll2007h , the mass of the gas is probably i ndeed 
negligible for their sample of galaxies. ICombes et Tall (2007) esti- 
mate the mass of the molecular gas for the same sample of galaxies 
and find masses of the order of some 10 7 Mq, w hich is about 3 to 
4 orde rs of magnitudes less than the results in Cappcllari ct al. 
||2006| ) suggest for the total masses of the galaxies. 



recall the metallicity gradient in elliptical galaxies, which is 
natural if they are more complex than MWGCs and thus 
more diverse in their internal properties. 



4 THE OBSERVED M/L v RATIOS AND 

PREDICTIONS OF STELLAR POPULATION 
MODELS 

For the remainder of this paper, we will compare the ob- 
served M/Lv ratios of the objects discussed in the previous 
sections to predictions from stellar population models, with 
the focus on the M/Lv ratios of the MCOs. 



4.1 The MCOs as simple stellar populations 

In order to find which stellar population models are appro- 
priate for the MCOs, we recall that most of the objects 
discussed here are old and note that a super-solar abun- 
dance of a-elements seems to be typical for the dynam- 
ically hot stellar sys tems discussed he r e, see e.g ICarnevi 
(1 19961) for MWGC s lEvstigneeva et ail (|2007l ) for MCOs, 
and Annibali et alj (|2007 ) for elliptical galaxies. Therefore, 
self-enrichment through the ejecta of type I supernovae 
(SNI) apparently does not play a major role in these sys- 
tems, as SNI are i mportant contributors of iro n to the in- 
terstellar medium (Ma tteucci fc Greggid Il986). This can 
be taken as an indicator for a stellar population with a 
narrow age spread, if the progenitors of SNI are assumed 
to be white dwarfs that surpass t he Chandrasekhar limi t 
by accretion of add i tiona l mat ter dWhelan fc: Ibenl Il973l ) . 
iMatteucci fc Recchil (|200ll ) and lGreggid l|2005l ) suggest me- 
dian time scales between some ten Myr and a few Gyr for 
the evolution of white dwarfs into SN I, depending on the 
initial conditions for the population. Considering stellar sys- 
tems with ages of ~ 10 Gyr, this can be taken as a rather 
short time scale. The assumption of populations of coeval 
stars within each stellar system thereby seems a reasonable 
approximation for at least MWGCs and MCOs. 

Besides age and age spread of the stars, a discussion of 
the M/Lv ratios of stellar systems has to account for the 
metallicities of their stars, since the metallicity is known to 
have a influence on the colour and the luminosity of a star 
with a given mass. Therefore the metallicities of the stars 
have to be known if one intends to construct a model for a 
stellar population which accurately describes a real stellar 
population, including its M/Lv ratio. 

In the following, two assumptions for the metal abun- 
dances in the stellar populations of the MCOs are made. 
This is not only for the sake of simplicity but also for the 
lack of more detailed data in most cases. 

Firstly, it is assumed that the metallicity-luminosity de- 
pendency of the stellar system can be characterised by the 
mean metallicity Z of the stellar system. This would cer- 
tainly be the case if Z was equal to the metallicities of the 
component stars, i.e. if all stars had the same metallicity. 
However, this is not necessarily the c ase for the stars in 
MCOs, as the examp les of u> Cen (e.g. iKavser et al. 1 120061 : 
IVillanova et alj|2007l ) and Gl (|Mevlan et aljhoOll ) show. On 
the other hand, imposing a more complicated metallicity dis- 
tribution on the stars of the unresolved stellar populations 
of the other MCOs does not seem reasonable. 



12 Dabringhausen et al. 



Table 3. MCOs with published metallicity estimates. [Fe/H] in 
Column 3 is taken as the measure for Z of the object. Also (V— 7) 
colour indices are given for some objects whose metallicities were 
derived from line indices. They provide the opportunity to test 
the validity of eq. ((HJ on a sample of MCOs (see section l5.21 . The 
columns of the table contain the following information: Column 
1: The name of the object, Column 2: [Z/H] if given in the ref- 
erence, Column 3: [Fe/H] either from the reference or calculated 
using eq. J7J, Column 4: The (V — I ) colour index, Co l umn 5 : 
T he reference to the sou rce paper: 1: lEvstigneeva et al] ll20 07t). 
2:lMieske et al.l |200rJ) 3: lHasegan et al. N2005I 1, 4: iMevlan et al] 
l|200ll l. 5: lHarrisl l|l996h . 



MCOs, we adopt [a/Fe] 
the relation 



0.3 for each one of them and use 



Name 


[Z/H] 




[Fe/H] 




(V-I) 


Ref. 


VUCD1 


-1.35.. 


. -0.33 


-1.12 


±0.51 


0.96 


1 


VUCD3 


0.00. . 


. 0.35 


-0.107± 0.175 


1.27 


1 


VUCD4 


-1.35. . 


. 0.33 


-1.12 


±0.51 


0.99 


1 


VUCD5 


-0.33. . 


. 0.00 


-0.447±0.165 


1.11 


1 


VUCD6 


-1.35. . 


. -0.33 


-1.12 


±0.51 


1.02 


1 


VUCD7 


-1.35. . 


. -0.33 


-1.12 


±0.51 


1.13 


1 


S417 


-1.35. . 


. 0.00 


-0.957 ±0.65 




1 


UCD1 






-0.38 


±0.05 


1.11 


2 


UCD2 






-0.90 


±0.33 


1.12 


2 


UCD3 






-0.52 


±0.11 


1.18 


2 


UCD4 






-0.85 


±0.29 


1.12 


2 


UCD5 














S314 






-0.50 






3 


S490 






0.18 






3 


S928 






-1.34 






3 


S999 






-1.38 






3 


H8005 






-1.27 






3 


Gl 






-0.95 


±0.09 




4 


lj Cen 






-1.62 






5 



Secondly, it is assumed that the mean iron abundance, 
[Fe/H], allows solid conclusions on Z. This assumption can 
be motivated with the finding that [a/Fe] ~ 0.3 seems not 
only to be true for MWGCs (|Carnevl li"996ri. but also for 
most of the MCOs that were analysed by Evstigneeva et all 
l|2007f ). This value appears to be very typical for massive, 
dense star clusters. 

The approximations and assumptions that have been 
made here and in Section 13.41 imply in their entirety that 
the stellar populations in MCOs can be considered as sim- 
ple stellar populations (SSPs), meaning that all stars and 
stellar remnants have the same age and the same chemical 
composition. 



4.2 The metallicities of the MCOs 

Informati on on the metallic it ies of MCO s are pub- 
lished in lHasegan et~al] (I2005TI. iMieske et al] I2 OO6). and 
lEvstigneeva et all ( 20071 ). lEvstigneeva et al.1 1 20071 ) give for 
each of the MCOs they examined an interval in which the 
actual mean metallicity Z of the MCO lies. We assume that 
this true val ue for Z of the MCO l ies in the middle of the in - 
terval given. lHasegan et all l|2005n and lMieske et al.1 (|2006l ) 
do not give estimates for Z of the objects they disc uss, but 
for [F e/H], Based on the obse r vation al findings by ICarnevI 
(1996) and lEvstigneeva et al] (|2007t ) and the assumption 
that the iron abundance characterises the metallicity of the 



[Z/H] = [Fe/H] + 0.94 [a/Fe] (7) 

found by iThomas et all (|2003T ) to calculate [Z/H] from 
[Fe/H]. The values that are adopted for the element abun- 
dances of the MCOs are summarised in Tab. [3] 

For the objects in Centaurus A no metallicities have 
been published so far, b ut (B — V) and ( V — I ) colour indices 
for them are available in iReikuba etall(|2007h . Observations 
show that there is a correlation between colour indices and 
[Fe/H] in GC systems. On this basis, an estimate of [Fe/H] in 
the objects in Centaurus A can be made by assuming that 
they follow a relation between colour and metallicity that 
has be en established on another GC system. iBarmbv et al] 
(2000) give relations between [Fe/H] and (V — I) as well as 
[Fe/H] and (B - V) for th e GC system of the Milky Way, 
using the data from lHarrisi (|l996l ): 

[Fe/H] (v _ /} = (4.22 ± 0.39) x (V - I) - (5.39 ± 0.35) (8) 
and 

[Fe/H] {s _ v) = (5.50 ± 0.33) x (B - V) - (5.26 ± 0.23). (9) 

The confidence range of these equations is set by the values 
(V — I) and [Fe/H] can assume for MWGCs. Their values 

for [Fe/H] are mostly between —2 and —0.5 dex. 

The advantage of the relations from IBarmbv et al] 

(2000) is that they have been established for both colour 
indices that have been measured for the objects in Centau- 
rus A, i.e. they allow us to fully benefit from the available 
data. Their disadvantage is that they do not account for 
a slight curvature in the relation between [Fe/H] and the 
colour indices, wh ich is typical for this relation according to 
lYoon et al] l|2006T ). However, given the apparent weakness 
of this departure from linearity, it seems justified to neglect 
it. 

We calculate [Fe/HLy_« and [Fe/HL s _ V j for each 
cluster in Centaurus A from eq. (JSj and (|9j if both colour 
indices are available. The results from eq. ([${ turn out to 
be systematically lower by sa 0.6 dex on average than the 
results from eq. ([9]), as can be seen in Fig. [6] It is obvious 
that the different results for the iron abundance calculated 
from different colour indices may indicate a serious problem 
with those estimates. A discussion on how reliable the re- 
sults based on these metallicity estimates are will be given 
Section [S] For now, we clearly distinguish between objects 
with [Fe/H] estimates from colour indices and objects with 
[Fe/H] estimates from line indices. 

The relation between [F e/H] and (V — I) colour found 
bv lKissler-Patig et al] (119981 ) by including (beside MWGCs) 
GCs around NGC 1399 has a slightly flatter slope than 
eq. (JHJ. It yields however similar results in the colour range 
interesting for the purpose here (deviations would be » 
0.2 dex in the most extreme cases). 

As a compromise between the two values that are es- 
timated for the iron abundances of the objects in Centau- 
rus A, we adopt the mean of both values as our final value 
for [Fe/H]. The error to this value has two components. The 
first of them is due to the intrinsic uncertainties to eqs. |8j| 
and @. The second component is the uncertainty due to the 
systematic difference between the results from eqs. (|8]) and 
JSJ. We estimate this error as half the difference between 
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Table 4. Colours and derived [Fe/ H] for the Centaurua A objects. The contents of the columns in the table are the following: Column 
1: Identification of the object like in lReikuba et aL | ll2007l) . Column 2: The (V - I) colour index, Column 3: The (B - V) colour index, 
Column 4: [Fe/H] calculated from the (V — I) colour index, Column 5: [Fe/H] calculated from the (B — V) colour index, Column 6: Our 
final estimate for [Fe/H] with the adopted errors. 



Name 


(B - V) 


(V - 1) 


[te/H\ (B _ v) 


l"nv, /til 
[te/HJ (v _ /) 


rev. /trl 


HGHH92-C7 


0.75 


0.91 


— 1.13 


— 1.55 


— 1 34 + 30 


HGHH92-C11 


0.94 


1.12 


—0.09 


—0.66 


—0 38 + 39 


HHH86-C15 


0.89 


1.03 


—0.36 


— 1.04 


—0 70 + 42 


HGHH92-C17 


0.77 


0.88 


-1.02 


-1.68 


-1.35 ±0.39 


HGHH92-C21 


0.78 


0.93 


-0.97 


-1.47 


-1.22 ±0.33 


HGHH92-C22 


0.79 


0.91 


-0.91 


-1.55 


-1.23 ±0.39 


HGHH92-C23 


0.76 


0.78 


-1.08 


-2.10 


-1.59 ±0.55 


HGHH92-C29 


0.89 


1.08 


-0.36 


-0.83 


-0.60 ±0.35 


HGHH92-C36 


0.73 


0.85 


-1.24 


-1.80 


-1.52 ±0.35 


HGHH92-C37 


0.84 


0.99 


-0.64 


-1.21 


-0.93 ±0.37 


HHH86-C38 


0.78 


0.91 


-0.97 


-1.55 


-1.26 ±0.36 


HGHH92-C41 


0.89 


1.09 


-0.36 


-0.79 


-0.58 ±0.33 


HGHH92-C44 


0.69 


0.85 


-1.47 


-1.80 


-1.63 ±0.26 


HCH99-2 


0.74 


0.84 


-1.19 


-1.85 


-1.52 ±0.39 


HCH99-15 




1.06 




-0.92 


-0.62 ±0.23 


HCH99-16 




0.79 




-2.06 


-1.76 ±0.23 


HCH99-18 


0.89 


0.89 


-0.36 


-1.63 


-1.00 ±0.67 


HCH99-21 




0.78 




-2.10 


-1.80 ±0.23 


R223 


0.80 


0.95 


-0.86 


-1.38 


-1.12 ±0.35 


R261 


0.83 


0.99 


-0.70 


-1.21 


-0.95 ±0.35 
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Figure 6. Comparison between [Fe/H]^_jj and [Fe/H]j fl _yj 
for the objects in Centaurus A. The numbers have been calcu- 
lated with eq. J8)l and eq. l[9} respectively. The dashed line indi- 
cates equality of [Fe/HL^^j and [Fe/H] ^ B _ v ^ . The dotted line, 
corresponding to [Fe/H]^_ J j = [Fe/H]^ g _ v ^ — 0.6, is a fit by 
eye to the actual distribution of the data. 



both estimates for a particular object. For the total error 
to the estimate of [Fe/H], the square root of the sum of the 
squares of both errors is assumed. 

For three objects only a (V — I) colour index is given. In 
these cases we simply set [Fe/H]^ v _ / j ±0.3 dex = [Fe/H], as 
0.3 dex is the average value by which the (V — /) colour 
indices of the other objects are changed. For estimating 



an error to these values for [Fe/H], the scatter of the 
data for [Fe/H], l/ _ / j and [Fe/H]^ fl _ V j around the relation 
[Fe/H]^ v / j = [Fe/H]^ B V ,j — 0.6 is calculated for the ob- 
jects in Fig [6] The scatter, s, is given by the equation 
* 2 = 7T=T £f [Pe/H] (v _ f)i - ([Fe/H] (s _ v) i - 0.6 dex)] 2 , 
where N = 17 is the number of objects in Fig. 6. This re- 
sults in s = 0.23 dex, which we adopt as the error to the 
[Fe/H] values of these three objects. 

The numbers for the metallicities of the objects in Cen- 
taurus A are listed in Tab. [4] 

Note that lHasegan et al.l (|2005T ) obtain the [Fe/H] es- 
timates for their objects also by comparison of the colour 
indices to the ones of GCs, i.e. in very much the same 
fashion as is done here for the objects in Centaurus A. 
The only MCOs with abundance estimates from line in- 
dices and thus estimates directly linked to an actual pres- 
ence of the according element s in the cluster ar e the objects 
from Evstigneeva et ail (|2007f l. the objects from lHilker et all 
|2007l ). uj Cen and Gl. 

Since Figs. [21 El H] and suggest a rather fluent transi- 
tion from the properties of MWGCs to the ones of MCOs, it 
seems worthwhile to include them in the discussion further 
on. A comprehensive co mpilati o n of t he iron abundances of 
M WGCs i s prov ided bv lHarrisl l|l996l ). Based on the results 
of ICarnevI (|l996t ). we assume [a/Fe] = 0.3 in order to calcu- 
late Z for them, as we did for the MCOs (eq. . 
Like the ones of MCOs, the stellar populations of MWGCs 
can be considered as old and coeval, but, in contrast to the 
ones of MCOs, dynamically evolved (i.e. loss of low-mass 
stars though evaporation driven by two-body relaxation). 

The [Fe/H] that are adopted for the MWGCs and the 
MCOs are plotted in Fig. [7] A tendency to higher abun- 
dances with higher masses is undeniable. Note however that 
selection effects might play a role here. There is a bias 
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Figure 7. The iron abundances adopted for this work plotted 
against the dynamical mass for the MWGCs and the MCOs. Open 
circles represent the MWGCs, filled circles the MCOs with abun- 
dance estima tes from line indices , open diamonds the MCOs in 
Centaurus A l lReikuba et al .12007ft and filled diamonds the MCOs 
in the Virgo cluster from lHaseean et al. The values for 

[Fe/H] of the latter two have been calculated from colour indices. 
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Figure 8. A comparison between the canonical IMF and the 
Chabrier IMF in the interval from 0.1 Mq to 10 M Q . Both IMFs 
are normalised such that J S,(m)mdm = 1, where m is the 
mass in solar units. The two IMFs are barely distinguishable on 
the whole mass interval. They actually are identical above a mass 
of lMg, except for a slightly different normalisation factor. 



against metal-rich objects for MWGCs, because they are 
concentrated towards the bulge of the Galaxy and th erefore 
harde r to observe than the metal-poor halo MWGCs (|Harris! 
1 19761 ). The GC systems of elliptical galaxies, on the other 
hand, have a larger fraction of red (probably metal-rich) 
GCs, which are also so mewhat brighter than the blue (prob- 
ably metal-poor) ones (jHarris et alj|2005 IWehner fc Harris! 
120071 ) . 



4.3 Predictions for M/Lv ratios from SSP models 

If information on the dependency of M/Lv ratio of a SSP 
on Z is combined with the estimates on the metallicity of 
the MCOs, it can be appraised what differences in TV are 
not due to differences in Z. Theoretic a l estim ates of Ty for 
different Z are taken from iMarastonl lj2005h for SSPs that 
form ed with a canonical IMF or a Salpeter-Massey IMF and 
from iBruzual fc Charlotl (|2003l ) for SSPs that formed with 
a Chabrier IMF. 

The canonical IMF is a continuous multi-power law, 



£k(to) 



(10) 



with ai = 1.3 for m < 0.5 M and a 2 = 2.3 for m > 0.5 Mq. 
It has been constrained after a decade-long study of various 
biases and found t o be consistent with all resolved stellar 
popu lations so far (|Kroupa et al.lll993l ; iKroupallioOll . 120021. 
l200Sl ). The Chabrier IMF is given for m < 1M as 



£c(m) oc — exp 
m 



(log(m/M Q ) -log0.08) 2 
0.9522 



(11) 



and equals the canonical IMF for m > 1 M up to a normali- 
satio n facto r. The transition at 1 Mq is continuous ijChabrierl 
l200ll . 120031 ). This IMF cannot be distinguished from the 
canonical IMF within the observational errors (Fig. [8} . To 



simplify matters, we will therefore also refer to the Chabrier 
IMF as the canonical IMF. The S alpeter-Massey IMF is 
a sin gle power law with a = 2.35 (|Salpeterlll955l : [Massevl 
1998). The SSP models used here have been obtained under 
the assumption that the IMFs are defined from 0.1 M to 
100 M©. 

Note that the upper ma s s lim it of the IMF as sug - 
gest ed by lWeidner fc Kroupal (|2004T ), lOev fc Clarke] (120051 ) 
and iFigerl (|2005l ) is higher than 100 M , but this does not 
have a mentionable affect on the expected M/L ratios of the 
SSPs discussed here due to the scarcity of high-mass stars 
in them. 

A lower mass limit of 0.1 Mq for the IMF neglects the 
existence of brown dwarfs. This is probably unproblematic , 
if one follows the argumentation bv lThies fc Kroupal l|2007f ). 
They suggest that star-like objects and brown dwarf-like 
objects are different populations and thus their frequencies 
ca nnot be descri bed by a single, continuous IMF as e.g. 
m iKroupal l|200ll) . The combined mass functions for brown 
dwarfs and stars which they find for star clusters in the 
Milky Way have many fewer brown dwarfs. Assuming a sim- 
ilar situation in the MCOs, brown dwarfs are not expected 
to contribute more than a few percent to t heir total mas s 
(opposed to m 10% for a mass function as in lKroupa|[200ll ). 

The ages that are c onsidered h e re for the SSPs are 9 Gyr 
and 13 Gyr. Note that IMarastonl |2005l ) distinguishes be- 
tween different horizontal branch morphologies, but this has 
a negligible impact on the dependency of Ty on Z of an old 
SSP. 

The benefit from us ing b oth the S S P mo dels from 
IBruzual fc Charlotl l|2003l ) and IMarastonl (|2005l ) although 
they cover the same ages and use (in principle) the same 
IMF is that different stellar evolutionary models have been 
used for constructing them. 

In order to make statements on the Ty of objects with 
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Table 5. Fit parameters for the interpolation formula for Ty 
to the data from th e SSP models. The SSP m odels are from: 1: 
iMarastonl ||2005| 1, 2: lBruzual fc Charlotl [120031) . 



Model 


a 


6 


c 


Ref. 


Salpeter IMF, 9 Gyr 


3.33 


0.82 


2.30 


1 


Salpeter IMF, 13 Gyr 


3.37 


1.20 


2.84 


1 


canonical IMF, 9 Gyr 


3.42 


0.42 


1.51 


1 


canonical IMF, 13 Gyr 


3.46 


0.79 


1.88 


1 


canonical IMF, 9 Gyr 


3.70 


0.23 


1.23 


2 


canonical IMF, 13 Gyr 


3.48 


0.55 


1.71 


2 



any Z, an interpolation formula that covers the whole Z- 
interval is needed. While it should be fai rly simple, it should 
also c losely fit the M/Ly ratios that iBruzual &: Charlotl 
|2003T ) and IMarastonl (|2005h find for specific metallicities. 
A function of the form 



F 4 ([Z/H]) = (a 



[Z/H] 



-b, ,Mq 

+ c )-j— » 



(12) 



where the index i distinguishes the different SSP models, 
fulfils these requirements well enough as Fig. [9]visualises. It 
can therefore safely be assumed that deviant estimates for 
TV are not due to an inadequate interpolation formula, but 
due to incorrect assumptions on the stellar population in the 
MCOs or to a failure of the SSP models. The parameters 
a, b and c found in least-squares fits are given in Tab. [5] 
Comparing these parameters for different SSP models with 
the canonical IMF reveals that they do not only depend on 
the assumed age o f the SSP, but also o n wh ethe r the SSP 
model s come from IBruzual fc Charlotl (|2003l ) or IMarastonl 
l|2005l ). This results in noticeably lower expectations for the 
M/Ly ratio from the SSP models from IBruzual fc Charlotl 
(2003), if co mpared to an in terms of age and IMF identical 
model from IMarastonl (|2005h . This proves the relevance of 
different stellar evolutionary models for the predictions from 
the SSP models. 

It should be mentioned that the value of Ty for the 
highest m etallicity was left out for the fit of eq. (|12[) to the 
data from IMarastonl |2005l ). because the omitted value was 
obtained by using a different stellar ev olution model than for 
the other data from IMarastonl (|2005t ). Moreover, excluding 
it results into a much closer fit of Fi([Z/~H]) to the remaining 
data, which already cover the metallicity range of the MCOs 
and the MWGCs. 

Note that stellar evolution only raises the M/Ly ratio 
of a stellar population. The M/Ly ratio of a 13 Gyr old SSP 
therefore provides an upper limit for the M/Ly ratio of a 
stellar population with a certain metallicity and IMF, since 
stellar populations cannot be much older according to the 
current estimates o n the age of the universe (13. 73^0^5 Gyr; 
ISpergel et ai1l2007h . 

If the stellar population of a star cluster with metallicity 
Z\ is similar to one of the modelled SSPs, one would expect 
Ty to be close to the prediction from eq. (|12[) for the M/Ly 
ratio at Z\: 



Fi\ Zl . 



If the stellar PDMF and the age of the star cluster is known 
(or assumed) to be similar to one of the SSP models that 
were introduced above and Z\ has been measured, Ty of a 
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Figure 9. The dependency of Ty on Z for different SSPs. The 
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part of each panel. The squares correspond to the models for 9 Gyr 
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old SSPs. The lines indicate the interpolation (eq. 1121 between 
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cluster which has the metallicity Z2, but is identical to the 
first one in all other respects can be estimated: 

T v \z 2 » - E rr J - x TvW- (13) 

The division by Fi\z 1 is imposed by the condition that 
the estimate for Ty must not be changed for Zi = Zi. 
~£y\z 1 /Fi\z 1 is the factor by which the theoretical predic- 
tion for the M/Lv ratio of a stellar system differs from the 
value that is observed. The multiplication of these numbers 
with Fi\z 2 is not necessary in principle, but it scales them 
by a constant M/Ly ratio such that the predicted M/Lv 
ratio from an SSP model with metallicity Z2 is expected to 
coincide with an observed value, if the model is appropriate. 

In order to eliminate the differences in Tv that are 
caused by differences in metallicity among the MCOs in the 
sample, we estimate Tv for them as it would be if they all 
had the same metallicity. This can be achieved by setting 
Z2 identical for all objects while using the measured Z for 
Zi in eq. (JT3J: 

T v ,n ^xTv, (14) 

where our (arbitrary) choice for Z2 is the solar metallicity, 
Zq. We refer to the M/Lv ratios calculated this way as the 
"normalised M/Lv ratios", Tv,n- Note that a comparison 
of a whole sample of values of observed M/Lv ratios to a 
single prediction for the M/Lv ratio of a SSP (as done in 
Fig. I10|l becomes possible that way. 

The values for Tv, n turn out to be quite insensitive to 
the actual choice out of the six sets of parameters a, b and 
c that encode different SSP models. This is due to the fact 
that the functions describing the dependency of Tv on Z are 
almost identical up to a scale factor for all the model pop- 
ulations that are considered here, i.e. the ratio Fi\z 2 /Fi\z 1 
is almost independent of the SSP model chosen. This means 
that the Tv,n that are calculated here are very likely to 
be good representations of the M/L ratios the MCOs and 
MWGCs would have if all their stars had solar composition, 
even if their PDMFs are different from all mass functions 
discussed here. 

However, the choice of the SSP model certainly has an 
impact on the prediction for the M / Lv ratio of a population 
that completely fulfils the assumptions made for the model: 
For different models, the predictions on such a population 
would be different by about a factor of Fi([Z/E]) /F 3 ([Z/E]). 

4.4 The normalised M/Lv ratios of the MCOs and 
the MWGCs 

The results for Tv,n of the MCOs and the MWGCs assum- 
ing different SSPs are presented in Fig. 1101 

The general distribution of the plotted points in all six 
panels of Fig. QXj] still closely resembles the distribution of 
those points in Fig. [5] which represent the same objects but 
with their observed M/Lv ratios. However, the increase of 
the M/L v ratios from the MWGCs to the MCOs is less 
pronounced once the effect of the metallicity on the lu- 
minosity has been accounted for, since the metallicities of 
the MCOs are usually somewhat higher than the ones of 
MWGCs (Fig. [7). 

There is a large spectrum of values for the Tv,n of 



the MCOs, ranging from « 2M©L~ 1 V to « ISMqL" 1 ^. 
However, most of them lie between ~ 3M0 Lg F and w 
7M© Lg 1 ^. This still covers a large range of values, but tak- 
ing into account that the Tv,n of individiual MCOs typi- 
cally also are uncertain within a range of ~ 2M0Lq 1 v to 
~ 4 Mq Lq v , it is not necessary to discuss physical rea- 
sons that could provide this scatter. However, two extreme 
outliers deserve more attention. 

The first one of them is the faint MWGC NGC 6535, 
which has Tv,n ~ 15 Mq Lq 1 ^ with large errors. As it is 
not only faint, but also fairly close to the galactic centre 
(the position is I = 27° 18', b = 10° 44' in Galactic coordi- 
nates), an accurate determination of its radius and velocity 
dispersion may be difficult due to the contamination with 
foreground stars. Moreover, its velocity dispersion has been 
derived from unpublished measurements. We therefore ex- 
clude it from Fig. [TOl 

The second outli er is the MCO S999 in the Virgo cluster 
(jHasegan et al.ll2005l ). which is the object with the largest 
Tv,n in all panels of Fig. 1101 If this rather high value 
is not due to a flawed measu rement, a scenario proposed 
by iFellhauer fc Kroupal (|2006h might offer an explanation. 
They proposed an enhancement of the M/ Lv ratio of MCOs 
by tidal interaction with the host galaxy. If this is indeed the 
case for S999, a faint envelope of stars may be detectable 
around it. It is noteworthy that this model can only provide 
an explanation for the Tv,n for a few MCOs out of a larger 
sample, as it requires quite specific orbital parameters. 

A comparison of the predictions of the SSP models 
with solar metallicity with the values for calculated Tv,n 
shows that the bulk of MWGCs and MCOs with masses 
< 2 x 10 6 Mq has lower Tv,n than it would be expected 
based on the assumed SSP models. Fig. [3] immediately re- 
veals that these star clusters have relaxation times well 
below a Hubble time, which means that they are dynam- 
ically evolved due to their age. This result is therefore 
in (at least qualitative) agre emen t with the predict ion by 
iBaumgardt fc Makinol (|2003h and iBorch et alj (|2Q07T ). who 
expect, based on their numerical simulations, the M/Lv ra- 
tio of a star cluster in a tidal field to be lowered by dynamical 
evolution for most of its lifetime (See Section f3.4p . 

The MCOs however have a strong tendency to higher 
M/Lv ratios compared to the theoretical prediction for a 
SSP with the canonical IMF, even for a 13 Gyr old popu- 
lation. There is only one SSP model, where in most of the 
cases the model expectation for Tv,n is higher than the ac- 
tual Tv,n of the massive MCOs. This is the model with a 
13 Gyr old stellar population which formed with a Salpeter- 
Massey IMF. For a 9 Gyr old population which formed with 
a Salpeter-Massey IMF, there seems to be agreement be- 
tween the model prediction for Tv and the actual Tv - How- 
ever, assuming that the IMF is truly universal and recalling 
that the stellar PDMFs of MCOs should still reflect their 
stellar IMFs as their dynamical evolution is slow, it can be 
concluded that the stellar population of the MCOs should be 
well described by a SSP formed with the canonical IMF. The 
Salpeter-Massey IMF deviates in the low-mass part strongly 
from the canonical IMF and can thus be ruled out if the 
above assumptions hold. 

It should be noted that the finding of observed M/Lv 
ratios being higher than the theoretical prediction from a 
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Figure 10. Normalised mass-to-light ratio, Tv,n, against mass for the MWGCs and MCOs based on the assumption that their stellar 
population can be described with SSP models. The origin of the SSP model, the assumed IMF and the assumed age of the SSP are 
given in the captions in each panel. The filled circles represent Ty n of MCOs with measured Z, the open diamonds represent objects 
for which Z was estimated from colour indices and open circles represent MWGCs. The dashed line indicates Fi\z & ; i- e - th e M/Ly ratio 
that the interpolation formula for the dependency of Yy on Z predicts for Zq, our reference metallicity. All points below that line have 
a lower Yy than the model predicts at their metallicity, all points above it exceed the model prediction. Natur ally Fi \z^ is very similar 
to the M/Ly ratios at Zq given in the actual models, where such a direct comparison is possible. llMarastonl d2005T) have data on Yy 
for Zq, Bruzual & Chariot (20030 use different grid points). 
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SSP model does not mean that the mass function of the 
chosen SSP model is inappropriate. Likewise, an agreement 
between the observed M/Lv ratios and the prediction from 
the SSP model does not mean that the assumed IMF is cor- 
rect. Consider for instance the presence of non-stellar black 
holes or non-baryonic DM in the MCOs, that lead to a rise 
of the M/Lv ratio unaccounted for by any SSP model. How- 
ever, in case the SSP model systematically overestimates the 
M/Lv ratios of a sample of clusters, the model is certainly 
not a good description for the stellar population of the clus- 
ters. 

Even if it is assumed that the MCOs only contain stars 
and stellar remnants, the significance of the tendency for 
higher Tv,n of the MCOs compared to SSPs whose IMFs 
agree with the canonical IMF should still be discussed. The 
case of a 13 Gy r old SSP with the canonical IMF from 
IMarastonl (|2005l ) is of special interest and will therefore be 
treated in detail, because this is the model where the devi- 
ation of the Tv,n calculated for the MCOs from the theo- 
retical expectation is the least pronounced. The values for 
Ty,n agree in fact with the prediction from the appropri- 
ate SSP model within the error for a large fraction of the 
MCOs, as can be seen in the middle right panel of Fig. 1101 
On the other hand, if taken as a sample, the MCOs which 
are more massive than 2 x 10 6 Mq still have a clear tendency 
for higher normalised M/Lv ratios than one would expect 
from the SSP model. 

A possibility to test whether a tendency is a significant 
deviation from an expecta tion is Pearson's test for the good - 
ness of fit, as it is found in lBhattacharvva fc Johnson! (ll977T l 
(see Appendix IA1|) . We apply this test on the MCOs more 
massive than 2 x 10 6 Mq under the assumption that their 
values for Tv.n would scatter just as much to higher values 
as to lower values compared to the prediction for Tv, n from 
an appropriate model. 

The result of the test is then that the probability for the 
found (or an even more one-sided) distribution of the values 
for Tv,n of the MCOs more massive than 2 x 10 6 Mq around 
the expected value for a 13 G yr old SSP with a canoni- 
cal IMF from IMarastonl (|2005r / is < 0.005. The hypothesis 
that this SSP model can fully describe the properties of the 
MCOs can therefore be excluded according to this test. 

The reliability of this result can be doubted, because it 
is not entirely clear whether the sample of the 31 objects, 
for which M > 2 x 10 6 M is fulfilled, is large enough to 
apply Pearson's test for the goodness of fit. Moreover, the 
objects with the more uncertain metallicity estimates from 
colour indices are included in this sample. 

We therefore also apply the sign test, as described in 
iBhattacharvva fc Johnsonl l|l977t ) (see Appendix lA2|) . on the 
13 MCOs with metallicity estimates from line indices. The 
hypothesis to be tested is that there is no significant differ- 
ence between their values for T^.n and the theoretical ex- 
pect ation assuming a 13 Gyr old SSP with a canonical IMF 
from IMarastonl (|2005l l. The probability that the Ty,n are 
larger than the theoretical expectation in 12 or more cases 
is 0.002 according to this test, i.e. it is highly improbable 
that the hypothesis is correct. 

Both statistical tests thus suggest that stellar popula- 
tion models cannot explain the M/Lv ratios as long as a 
canonical IMF is assumed, even for the maximum age the 
stellar population could have in order to be consistent with 




.-| 2 1 1 1 1 1 1 1 1 

10 2,5 10 3 10 3,5 10 4 10 4 - 5 10 5 10 5 - 5 10 6 

mass [ 10 6 solar units] 

Figure 11. The 16th percentiles (lower open circles), the me- 
dian values (filled circles) and the 84th percentiles (upper open 
circles) of the distributions of the metallicities of galaxies in differ- 
ent total-stellar- mass bins (Gallazzi ct al. 2005). The lower dot- 
ted line, the solid line and the upper dotted line are our fits of 
eq. H15I I to the 16th percentiles, the median values and the 84th 
percentiles, respectively. 



Table 6. Best-fitting parameters of cq. H15I I if fitted to the me- 
dians (50th percentiles) of the distributions of the metallicities of 
galaxies in different mass bins, as well as to the 16th and 84th 
percentiles of these distributions . The required data o n the metal- 
licity distributions is taken from Gallazzi et al. (2005), their table 
2. 



Percentile 


a 


b 


c 


d 


Median (P50) 


0.29 


3.06 


-4.09 


-0.267 


P16 


0.41 


2.72 


-4.37 


-0.555 


P84 


0.13 


1.78 


-4.07 


0.118 



the age o f the universe accor ding to cosmological models. 
Note that Mie ske et al.l l|2006h suggest intermediate ages for 
the MCOs in the Fornax cluster. The actual discrepancy 
between the true values for M/Lv ratios and the SSP mod- 
els with the canonical IMF would then be larger than in 
the case discussed above. This means that as long as the 
SSP models do not fail to describe real stellar populations, 
the MCOs either contain additional non-luminous matter, or 
their PDMFs must be different from what one would expect 
for a stellar system formed with the canonical IMF. 

4.5 The normalised M/Lv ratios of elliptical 
galaxies 

We now compare the M/Lv ratio of elliptical galaxies and 
galactic bulges with the prediction for the M / Lv ratio of 
a 13 Gyr old SSP with the ca nonical IMF according to the 
models from IMarastonl |2005h . 

The metallicity estimate that enters the calculation of 
the normalised M/Lv ratio of the elliptical galaxies and 
galactic bulges is based on results on the metallicities of 
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Figure 12. The normalised M/Ly ratios of all objects plotted in Fi g. \Ej The s ymbo ls are as in Fig. 2. The stellar population model 
assumed for calculating the normalised M/Ly ratios is the one from iMarastonl j2005t ) for a 13 Gyr old SSP (i.e. the same SSP model 
as for the middle right panel of Fig. HOP . Black bars indicate for five of the elliptical galaxies the range of normalised M/Ly ratios they 
would assume if their metallicity would vary between the adopted values for the 16th percentile and the 84th percentile of the distribution 
of the metallicities of galaxies with that mass. The dashed line corresponds to the normalised M/Ly ratio expected according to the 
SSP model assumed for the objects in this figure. 



galaxies from the S loan Digital Sky Survey obtained by 
IGallazzi et all (|2005h . It is apparent from their data that 
the metallicities of galaxies in a given total-stellar-mass bin 
are distributed over a range of possible values (their figure 8 
and table 2). In the present paper, the median of this distri- 
bution is taken as a representative value for the metallicities 
of the galaxies in that mass bin. The metallicities of the ellip- 
tical galaxies and galactic bulges in our sample as a function 
of their mass are calculated using the function 



[Z/B](M) = a arctan b 



lop 



M 



10 6 M Q 



+ c 



+ d (15) 



with parameters a, b, c and d found in a least-squares 
fit to the median metallicities of galaxies in total-stellar- 
mass bins between ~ 10 9 M Q and « 10 12 Mp), as given by 
IGallazzi etall (|2005h . The data from IGallazzi et all (|2005T ) 
as well as the fit to them is shown in Fig. 1111 The best-fitting 
parameters a, b, c and d are noted in Tab. H3 

For the abundances of dSphs, it is assumed that 
their values for [Fe/H] can be identified with their val- 
ues for [Z/H]. Iron ab undance s for most dSphs dis- 
jMateol l| 19981), except for And 



cussed here are given in 



II dMcConnachie et alll2005l'). And XI |McConnachie et all 



2005) and UMa I jSimon fc Gehall2007T l 



The normalised M/L ratios which are implied by the 
adopted metallicities for the elliptical galaxies, the galac- 
tic bulges and dSphs introduced in Section [2] are plotted 
together with the normalised M/Ly ratios of MCOs and 
MWGCs in Fig. E] 

IGallazzi et all (|2005T ) find especially for low-mass galax- 
ies a large spread for the distribution of their metallicities. 
To quantify the uncertainties that arise for the adopted nor- 
malised M/Ly ratios from the range of likely actual metal- 
licities of galaxies, eq . (|15|l is also fitted to the values from 
IGallazzi et all (|2005h for the 16th and 84th percentiles of 
the distributions of metallicities of galaxies in different mass 
bins. The best fitting parameters a, b, c and d can be found 
in Tab. [6] Using these parameters, likely values for a high 
and a low metallicity in a given galaxy can be estimated 
depending on its mass and the according normalised M/Ly 
ratio can then be calculated. In Fig. 1121 the possible range 
of normalised M/Ly ratios suggested by the lower and the 
upper estimate of its metallicity is indicated for five sample 
objects with black bars. 

It thereby becomes apparent in Fig. [12] that the spread 
of the normalised M/Ly ratios of elliptical galaxies and 
galactic bulges cannot be explained by different metallicities 
alone, but that at least one more parameter (e.g. the mean 
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age of their stellar populations) must vary among them as 
well. 

Consider the elliptical galaxies and galactic bulges with 
the highest normalised M/Ly ratios. Given the adopted 
range for their likely metallicities, the range of M/Ly ra- 
tios possible for them is inconsistent with the prediction for 
their nor malised M/Ly rat io from a model for a 13 Gyr old 
SSP from IMarastonl (|2005h : especially for the objects with 
high dynamical masses. This suggests, as for the MCOs, 
an IMF different from the canonical IMF for their stellar 
populations or the presence of additional (gaseous or non- 
baryonic) matter in them. 



5 DISCUSSION 

5.1 How reliable are the SSP models? 

The results that have been obtained in Section 14.11 are 
strongly based on the reliability of SSP models which are 
in turn based on the reliability of evolutionary stellar mod- 
els. However, the reliability of these models cannot be taken 
for granted, as th e diffe rence s between the m odels from 
iBruzual fc Charlotl (|2003r i and IMarastonl (120051 ) already in- 
dicate. 

Another issue that may hint at difficulties with the SSP 
models is the relation between the iron abundance and the 
colour indices they suggest. This becomes apparent when 
using them to predict [Fe/H] of the MCOs in Centaurus A 
from their colours. This can be done by setting up alter- 
native equations to eqs. (JSj and ((9]l by fitting interpolation 
functions to the (V — J)-[Fe/H] value pairs and the (B — V)- 
[Fe/H] value pairs given by the SSP models (i.e. as in Section 
I4.3l for a relation between the metallicity and the M/Ly ra- 
tio). Fig. H3l shows that a good fit between the data and the 
interpolation can be achieved with functions of the form 

[Fe/H] SS p = a (V - 1) + b(V - /) ' 5 + c (16) 

for the (V — I) colour index and analogous for the (B — V) 
colour index. The subscribt SSP in eq. (|16|l is supposed to 
indicate that these estimates for [Fe/H] from colour indices 
are based on SSP models, in contrast to the estimates for 
[Fe/H] from eqs. © and (0, which are based on observations 
of the MWGCs. 

In Fig. 1141 [Fe/H] (v-/), ssp is plotted against 
[Fe/H]( S _v),ssp f° r the objects in Centaurus A. Each 
panel represents a choice of the SSP model which is 
assumed to represent the stellar population of the objects 
in Centaurus A best. There are two features of the dis- 
tribution of the data, which are remarkably little affected 
by that choice. The first one is the undeniable tendency 
for [Fe/H] (V-/),ssp < [Fe/H]( S _v),ssp- The second one is 
that the spread of the values for [Fe/HW_n ggp is larger 
than the spread of the values for [Fc/H]( S _v),ssp- However, 
if one of the SSP models is an adequate description for 
the actual SSPs in Centaurus A, no systematic difference 
between the two estimates for [Fe/H] from this SSP model 
would be expected. 

One could therefore come to the conclusion that none 
of the SSP models considered in this paper reflects the ac- 
tual stellar populations of the objects in Centaurus A. Note 
however that neither assuming an age of 5 Gyr nor consider- 
ing a different horizontal branch morphology for the models 



from lMara ston (2005]) can enhance the concordance between 
[Fe/H]( fl _v),ssp an d [Fe/H](v-j),ssp for the objects in Cen- 
taurus A. This could be evidence of the standard SSP mod- 
els failing to give a detailed an d accurate descrip tion of real 
stellar populations in principle. IXin et all (12007! ) claim that 
this might indeed be the case as long as SSP models are only 
based on the evolution of single stars but neglect the exis- 
tence of blue stragglers, which are thought to be products of 
stellar interactions. Given the complex abundance patterns 
in resolved massive star clusters, it also seems well possible 
that the observed (integrated) (B — V) and (V — I) color in- 
dices of the objects in Centaurus A can only be reproduced 
by stellar population models which account for an age and 
metallicity spread of the stars. 

An alternative explanation for the inconsistency be- 
tween [Fe/H]( S _^),ssp an d [Fe/HW_n gsp could be a so far 
unidentified observational bias in the colour observations of 
the objects in Centaurus A. This notion is made attractive 
by the finding that applying the observed relations eqs. (JSj 
and ((9J for the estimation of [Fe/H] leads to smaller [Fe/H] 
estimates from (V — I) colour indices than from (B — V) 
colour indices for the MCOs in Centaurus A as well (Fig. 
[6]). If the difference between the metallicity estimates from 
eqs. © and was, for instance, caused by a systematic 
error to the (B — V) colour indices, their offset from the 
true (B — V) colour indices would be « 0.1 dex. 

Considering both the inconsistency of the iron abun- 
dances derived from the different colour indices by using the 
SSP models and the noticably different predictions of differ- 
ent SSP models on the M/ Ly ratio of the same population, 
it still seems possible that the enhancement of the M/L ra- 
tios of the MCOs compared to the theoretical predictions for 
SSPs with the canonical IMF is due to a failure of the SSP 
models. 

5.2 How reliable is an estimate of [Fe/H] from 
colour based on observations? 

The alternative to estimating [Fe/H] from colour indices 
based on a SSP model is the approach chosen for this paper, 
namely using a relation between [Fe/H] and colour indices 
that has been established on a sample of observed star clus- 
ters, such as eqs. © and ©. But just like the estimate of 
[Fe/H] by using SSP models, this approach is not unprob- 
lematic, as will be discussed here. 

It is helpful to define two terms for the further discus- 
sion: We call the sample of objects for which the relation 
between [Fe/H] and colour was established the "calibration 
sample" . The sample for which only colour indices are mea- 
sured and where the relation between [Fe/H] and colour is 
used for a metallicity estimate is called the "target sample" . 
In our specific case, the MCOs in Centaurus A are the tar- 
get sample and applying eqs. (|8]) and <(9j on them makes the 
MWGCs the calibration sample. 

There are two problems, that are generally attached to 
an estimate of the iron abundances from the colours of ob- 
jects in a target sample based on observations of an calibra- 
tion sample. Firstly, it has to be assumed that the objects 
in both samples have at least typically the same PDMFs 
for shining stars and the same ages. If this is not the case, 
this method is likely to fail because colours depend on these 
parameters as well as on metallicity. 
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Figure 13. The relation between colour indices and [Fe/H] according to SSP models. The right panels show [Fe/H] against the (V — I) 
colour index while the left panels show [Fe/H] against the (B — V) colour index. Squares show the data for 9 Gyr old populations. The 
dashed line is the fit to them. Circles show the data for 13 Gyr o ld populations. The d otted line is the fit to the m. The thin solid lines 
represent the relations that have been established for MWGCs bv lBarmbv et al.1 goOO). Out of the SSP models bv lMarastonl J2005T) . the 
case of a red horizontal branch is shown in this figure. This morphology is said to reflect the horizontal branches in most of the metal-rich 
GCs and therefore seems to be an appropriate choice for the MCOs, which show similar metallicities if measured. 



Secondly, relations such as eqs. © and © are only 
fitting formulae to a data sample with scatter. However, if 
these relations are applied to the objects in the calibration 
sample, the resulting estimates for [Fe/H] lie in the same 
parameter space as the values for [Fe/H] from line indices. 
The same is true if the calibration sample and the target 
sample are indeed comparable. 

As a test whether the MWGCs are a good choice for the 
calibration sample for the MCOs in Centaurus A, the values 



for [Fe/H]^ v _ 7 j from eq. (]5J) are compared to the values for 
the estimates of the iron abundances from line indices (and 
thus directly linked to a observed iron content in the star 
clusters), [Fe/H] obs . The published data (see Tab. [3} allow 
such a comparison for the ten objects plotted in Fig. 1151 

There is no significant trend for [Fe/H]^ v _ J j to be larger 
or smaller than [Fe/H]„ Kp , as the application of the sign test 
( Bhattach arvva fc Jo hnson 1977;, Appendix lA2|) shows. Un- 
der the hypothesis that there is no significant difference be- 
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Figure 14. A comparison between [Fe/H]( B _^)gsp and [Fe/H] i"),SSP f° r the objects in Centaurus A (Tab. [2] and Tab. [4}. The 
values are estimated by using the fits to the data from the SSP models plotted in Fig. 1131 The dashed line indicates equality of both 
estimates for the iron abundances. If there was no systematic difference between them for the objects in Centaurus A, the distribution 
of the data would follow these lines. Errors to the plotted points arc not shown. They are probably governed by the errors to the colour 
measurements (which are unknown to us) and by a mismatch between the SSP models and the real stellar population s of MCOs (whic h 
is to be shown by this figure), but not by the errors to the interpolations plotted in Fig. 1131 For the SSP models from iMarastorj (2005), 
a red horizontal branch is assumed. However, this does not have a strong impact on the distribution of the data in the according panels 
of this figure. 



tween the two values, the probability for having only four 
or less out of ten with [Fe/H] obs > [Fe/H] (y _ I) is 0.377. A 
result as the one plotted in Fig. [15] is therefore quite prob- 
able. From this point of view it seems justifiable to apply 
eq. iJHJ on the MCOs, although it was originally fitted to the 
MWGCs. 

Recall however that [Fe/H], s v j is systematically 



higher than [Fe/H]^_ 7 j for the objects in Centaurus A 
(Fig. HJ}. Since we adopt the mean of [Fe/H]^ „ and 
[Fe/H]^ fl _ v ,j, [Fe/H] of the clusters in Centaurus A will be 
overestimated if [Fe/H] ( - v _ J j reflects their true abundances 
well. This is a conservative choice in our case, because a 
higher estimate for [Fe/H] leads to a lower estimate for Ty,„. 
We arrived at the result that a SSP model with the canonical 
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Figure 15. [Fe/H]/y_j\ calculated from eq. [8] plotted against 
the iron abundance estimate from line indices, [Fe/H] obs , for the 
objects whose (V — I) colours are given in Tab. \3\ The dashed 
line indicates equality between observed and calculated values. 
Apparently there is no clear tendency for the points to be lo- 
cated only on one side of the line. This indicates that there is no 
systematic difference between [Fe/H](y_j-j and [Fe/H] oba . 



IMF underpredicts the M/Lv ratios of the MCOs neverthe- 
less, this therefore being a robust conclusion. 



5.3 The impact of a wrong estimate of [Fe/H] on 
the comparison of the dynamical M/Lv ratios 
with the SSP models 

As the metallicities of the MCOs may be subject to system- 
atic errors, it makes sense to discuss the impact of a wrong 
metallicity estimate on our claim that the M/Lv ratios of 
the MCOs are inconsistent with the predictions from SSP 
models with the canonical IMF. We discuss one case in detail 
in order to give an impression how this affects our results. 

Suppose the objects in the calibration sample are well 
described by a SSP with the same mass function, but that 
the target sample is younger than the calibration sample. 
The colour of the objects in the target sample is then bluer 
than it would be if they were of the same age as the objects 
in the calibration sample. 

When relations like eqs. (JSJ) and @ are applied in order 
to estimate the iron abundance, it is implicitly assumed that 
the stellar populations of the objects in the target sample 
are the same as the ones in the calibration sample. The 
estimates for the iron abundances are therefore too low if 
the target sample is younger than the cal ibration sample , 
because of the age- metallicity degeneracy (|Worthevl ll994). 
As a consequence, the Tv,n calculated from eq. (fl4|) is too 
high, since the denominator on the right side of eq. (|14|l only 
decreases with decreasing Z due to the exponential nature 
of eq. ((12) . 

However, the prediction for Tv,n made by a SSP model 
increases with the assumed age of the SSP. The expectation 
for the Tv, n of the objects in the target sample is therefore 
also too high, if they are compared to an SSP model which is, 



concerning the assumed age of the objects, more appropriate 
for the objects in the calibration sample. Thus, the error 
that is made in the estimation of the values for [Z /~H] of the 
objects in the target sample by assuming a common age for 
all objects tends to balance the error that is made when all 
objects are compared to the same SSP model. 

An analogous argument can be found if the objects in 
the target sample are depleted in low-mass stars compared 
to the objects in the calibration sample. In this case, it is 
the scarcity of low-mass stars that makes the objects in the 
target sample bluer and thereby leads to a too low metal- 
licity estimate for them. The resulting too high estimate for 
Tv'.n for these objects is compensated if they are compared 
to a SSP model with a full population of low-mass stars 
(which are faint and therefore enhance the M/Lv ratio of 
the stellar population). 

The reverse argumentation can be applied to objects 
with a higher age or more low-mass stars than the objects 
in the calibration sample. 

It thereby seems that, also for objects with metallicity 
estimates from colour indices, finding the values for Tv,n 
above the expectation from the SSP model really is an indi- 
cator for additional non- luminous matter in the object. 



5.4 Implications of a high M/Lv-r&tio in the 
MCOs 

Two explanations for the systematic enhancement of the 
M/Lv ratios of the MCOs more massive than 2 x 10 6 Mq 
compared to the predictions from SSP models with the 
canonical IMF are possible. 

The first possibility is that the mass ive MCOs are 
embed ded in DM haloes, as proposed by lHasegan et al.l 
(2005). However, for MCOs with small effective radii and 
high M/Lv ratios, the mean density of the DM within five 
half-light radii would have to be between 1 Mq pc -3 and 
10 Mq pc -3 in order to have the observed impact on their dy- 
namics. Adopting the universal DM density profiles as they 
are p redicted by standard A CDM cosmology dNavarro et al.l 
1996), only DM haloes with masses of 10 12 M or more could 
accumulate enough DM in their central parts (Dabring- 
hausen & Kroupa 2008a, in preparation). 

The alternative to suggesting non-baryonic DM in the 
MCOs is to give up the notion of a universal IMF for all 
stellar populations. Such an alternative IMF would either 
be over-abundant in low-mas s stars with high M/Ly ra- 
tios (bottom-heavy IMF), see iMieske fc Kroupal (|2008h . or 
over-abundant in massive stars (top-heavy IMF). The lat- 
ter possibility would imply a high number of dark stel- 
lar remnants in an old stellar population. Especially a 
top-heavy IMF seems attractive, since it is also suggested 
by models for galaxy evolution (e.g. [ Baugh et al.l 120051 ; 
iNagashima et al.1 [ioOol; Ivan Dokkurri 120071) or GC evolu- 



tion (e.g. iD'Antona fc Caloil |2004| ; iPrantzos fc Charbonnell 
2006). These issues will be examined in more detail in 
a forthcoming paper (Dabringhausen & Kroupa 2008b, in 
preparation). 

5.5 On the nature of MCOs 

Apart from the finding that the MCOs more massive than 
2 x 10 6 Mq are in disagreement with the expectations for 
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their M/Lv ratios according to SSP models with the canoni- 
cal IMF, the increase of typical radii at about the same mass 
is probably the most intriguing observation. This raises the 
question whether the massive MCOs (mostly classified as 
UCDs in the literature) constitute a population different to 
other populations of stellar systems as far as their origin is 
concerned. This question is of special interest for the rela- 
tion between massive MCOs and GCs, since the seemingly 
continuous rise of the mean radius above 10 6 Mq makes the 
notion of a single population of objects attractive (single 
population in the sense of a common scenario that leads to 
their formation). In this case, the evolution of such an ob- 
ject must be different at very high masses in order to account 
for the increase of radius with mass. A possible reason for 
this could be a dependency of star formation on gas density 
(which increases with mass for objects with the same exten- 
sion, consider e.g. the MWGCs in Fig. [4j) . If this dependency 
would lead to a greater mass loss during the lifetime of the 
cluster, it could explain the greater extensions because mass 
loss enlarges a star cluster. A possible example would be a 
top-heavy IMF in very dense star forming regions, which 
would cause a stronger mass loss by type II supernovae. 

On the other hand, there have been efforts to de- 
sign models that can specifically reproduce the parame- 
ters of UCDs. One of th ese models is the scenario by 
iFellhauer fc Kroupal (|2002h that UCDs are the merger of 
massive cluster complexes as are seen to be forming in mas- 
sively interacting galax ies. Another o ne is the scenario by 
iBekki et al] {200$ ) and lGoerdt et al.1 (|2008h that UCDs are 
the cores of nucleated galaxiefl Under the condition that 
GCs form in the collapse of a single molecular cloud, ob- 
jects that were formed in one of the above scenarios would 
indeed be of a different origin. This would offer natural ex- 
planations for the masses and the radii of those objects to be 
larger than for typical GCs. In this case, UCDs and GCs are 
two different populations that mix in the mass interval from 
10 6 M© to 10 7 Mq because both kinds of objects formed in 
intense starbursts that converted a similar amount of gas 
into stars. 

A question connected to the issues discussed here is 
whether it is expedient to discriminate the MCOs into UCDs 
and GCs. We think that this distinction can be justified. It 
clearly makes sense if UCDs really formed in a different way 
than GCs. But it also makes sense in the case that GCs and 
UCDs were initially formed in the same way, thus are in 
principle to be considered as the same class of objects as far 
as their origin is concerned. In this case, "UCD"would be a 
useful term to emphasise the peculiarities, for example the 
higher relaxation times, which very massive clusters usually 
show in comparison with their low-mass counterparts. Thus, 
UCDs could be defined as those compact stellar systems, 
which have relaxation times longer than a Hubble time and 
thereby are (almost) collisionless systems on this ti me scale. 
This d efinition is the same as the one proposed by iKroupal 
(1998) for a distinction between star clusters and galaxies, 
i.e. UCDs are galaxies in that sense. 

Also Gl and ui Cen are classified as MCOs instead of 



5 Note that the scenario Bckki et al. I <2003ft proposes is inconsis- 
tent with A CDM theory, because it has to assume that the DM 
haloes of the progenitors of the UCDs are cored instead of cusped. 



GCs in this paper, because of the spread that their stars 
show in [Fe/H] (w Cen certainly and Gl presumably). Pecu- 
liarities in element abundances can in principle be another 
way to discriminate UCDs from GCs by observational pa- 
rameters. However, the MCOs in other galaxy clusters can- 
not be resolved into stars with the current instrumentation. 
A similar pattern of the chemical composition of their stars 
as in the MCOs in the Local Group can for this reason only 
be presumed so far, but not be proven in the near future. 



6 CONCLUSIONS 

In this paper, a sample of compact stellar systems covering 
the transition from globular clusters (GCs) to ultra compact 
dwarf galaxies (UCDs) and referred to as massive compact 
objects (MCOs) in this work, is compared to other dynami- 
cally hot stellar systems. Moreover, the M/Lv ratios of the 
MCOs and the Milky Way GCs are compared to predictions 
from models for stellar populations. Our main conclusions 
are as follows. 

Departing from radii typical for GCs, which are con- 
stant with mass, greater extensions are correlated with 
higher masses for dense stellar systems more massive than 
10 6 M Q . A strong increase of the median two-body relax- 
ation time with mass is the natural consequence. We also 
find that stellar densities peak at a mass near 10 6 M©. 

Dwarf spheroidal galaxies (dSphs) take on a special po- 
sition among dynamically hot stellar systems. This is espe- 
cially apparent from their dynamical M/Lv ratios, which 
are in some cases higher by one to two orders of magnitude 
than for any other dynamically hot stellar system. Also note 
the large spread of the M/Lv ratios of the dSphs, which 
would imply very different DM densities in the visible parts 
of different dSphs, if the dSphs were in dynamical equi- 
librium. It therefore seems improbable that the masses of 
dSphs can be determined by simple application of Jeans' 
equations. 

The fact that compact stellar systems with t re i < th 
mostly have a much lower M/L ratio than systems with 
trei > th is qualitatively c onsis tent with th e findi ngs by 
iBaumgardt fc Makinol (|2003T ) and lBorch et~ai1 (|2007l ). How- 
ever the change to the M/L ratio by relaxation-driven mass 
loss does not exceed 0.3 Mq Lq 1 for most GCs, which is too 
little to explain the observed differences to the M/L ratio of 
GCs and UCDs. Dynamical evolution is slow for UCDs, as 
their high relaxation times indicate, and consequently the 
decrease of the M/Lv ratio by this process is slow as well. 
Moreover, the slow dynamical evolution leads to the stel- 
lar present-day mass function being almost identical with 
the stellar initial mass function for main sequence stars. We 
also found that the assumption of a population of old coeval 
stars in each massive MCO probably constitutes a good ap- 
proximation to their real stellar populations. 

Taken together, the lack of dynamical evolution and the 
narrow age spread of the stellar populations make a com- 
parison between the MCOs and theoretical predictions from 
SSP models with widely used IMFs reasonable. The SSP 
models also allow to account for the differences due to the 
different metallicities of the MCOs. The limiting factor here 
is, if the reliability of the SSP models is taken for granted, 
the only rough knowledge of the element abundances in the 
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MCOs. It turns out that the dynamical M/Lv ratios of the 
MCOs more massive than 2 x 10 6 Mq have a significant ten- 
dency to be even higher than the predictions of models for 
very old stellar populations, provided the IMF is chosen in 
agreement to the observations of stellar populations, where 
at present times low-mass main-sequence stars can be re- 
solved (i.e. populations in the Milky Way and in objects in 
its immediate surroundings, such as the Magellanic Clouds). 

It was shown however, that the SSP models that were 
used for the estimate of the expected M/Ly ratio of the 
MCOs cannot produce consistent [Fe/H] estimates for the 
objects in Centaurus A from the different colour indices mea- 
sured for them. This poses the question whether the SSP 
models in their current state (e.g. without binary evolution) 
are truly reliable. On the other hand, if the predictions for 
the M/Ly ratios from the SSP models are correct, the dis- 
crepancy between them and the dynamical M/Ly ratios ob- 
served in the MCOs suggests that the more massive MCOs 
contain DM or that the stellar IMF in some stellar systems 
is different to the ones of resolved stellar populations. Both 
possibilities will be studied in follow-up papers. 

Summarising, « 10 6 Mq is a critical mass-scale at which 
the system length-scale begins to increase, the highest stellar 
density is reached, the relaxation time becomes comparable 
to a Hubble time and evidence for dark matter appears. 
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APPENDIX A: STATISTICAL TESTS 

Al Pearson's test for the goodness of fit 

Pearson's test for the goodness of fit 
( Bhat tacharvva fc Johnson] ll977T ) can be used for de- 
ciding whether the frequency of a certain result for a 
measurement that has been performed on n objects de- 
viates significantly from an expected frequency. For the 
special case that only two results A and B can be the 
outcome of each measurement (A could be for example 
a result higher than a theoretical expectation and B the 
opposite case), result A will have occurred j times and 
result B n — j times. The probability of this outcome can 
now be calculated if a certain probability p for the case A as 
the result of a measurement is assumed. A useful measure 
for this is given by the equation 



ti-pn) 2 , ((n-j) - (1 -p)n) 2 



pn 



(1 -p)n 



(Al) 



There are tabulated values for \ 2 (e.g. ta ble 6 in the ap- 
pendix of [Bhattacharvva fc Johnsonl 1 19771) which make it 
possible to read off the probability for x 2 being higher than 
some value for a series of measurements, if the hypothesis 
for the probability p is correct. (The degree of freedom is 
one in this case.) 



A2 The sign test 



The sign test (iBhattacharvva fc Johnsonl 1 1977h is specif- 
ically designed for a small number of pairs of values, 
( X\ , X2 ) , and is supposed to detect whether there is a sig- 
nificant trend for X2 being larger or smaller than X\ or 
not. X\ and X2 could be two measurements under different 
conditions (e.g. other instruments), or X\ could be a value 
inferred from an observation, while X2 is the theoretical pre- 
diction for this value. If the conditions under which X\ and 
X2 were obtained do not result into systematically larger or 
smaller values for X2 compared to X\ , the probability for X\ 
being larger than X2 is 0.5. The probability that X2 > X\ 
for j out of n pairs of values is then given by the binomial 
distribution. 



